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S1!MB0LS 


Lower Case Latin 

a semi -diameter perpendicular to Z-axis in ellipsoid 

b semi-diameter parallel to Z-axis in ellipsoid 
f distributed loads in local coordinates 
h shell thickness; face sheet thickness 

i index; beginning edge of shell segment; independent Joint of 
kinematic link; subscript "inside" 

J index: ending edge of shell segment; dependent Joint of 

kinematic link 

n index on harmonic 

o subscript "outside" 

r radius 

s index of segment; coordinate in cylinder or cone 

t core thickness in sandwich shell 

w normal deflection, positive inward 

Upper Case Latin 


C stiffness eccentricity parameters; offset distance in ogive, ellipse 
D bending stiffness parameters 

E Young’s modulus (ib/in ) 

F lineal force (ib/in) 

G shear modulus (ib/in ) 

K extensional stiffness parameters 

M bending moment on shell (in-lb/in) 

N membrane stress resultant (ib/in) 

Q transverse shear stress resultant (ib/in) 

R radius; "global" coordinate, positive radially outward 
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SYMBOLS (continued) 


T temperature j "global" coordinate, tangential 

X Cartesian coordinate, 0 = 0 at X-axis 

Y Cartesian coordinate 

Z Cartesian and "global" coordinate, coincides with axis of revolution 

Greek 

a angle between rotated coordinates 

P ratio of semi -diameter parallel to Z-axis in ellipsoid to semi- 

diameter perpendicular to Z-axis 

7 shear strain; non-linear parameter; angle of inclination of 
kinematic link 

? normal coordinate, positive inward 

0 circumferential angular coordinate (rad) 

1 shell parameter 

V Poisson's ratio 

a normal stress (ib/in ) 

2 

T shear stress (ib/in ) 

meridional angulai* coordinate (rad) 
w rotational displacement (rad) 

Q rotational displacement in ^global” coordinates (rad) 

A displacements in fixed or "global” coordinates 

A segment length parameter 

Mi scellaneo us 

eq equivalent 

s4> sin <t> 

c<p cos 4) 

Other symbols are defined in the text where used^ 
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SECTION 1 


PROGRAM CAPABILEK 


The use of an accurate shell theory to analyze structural shell prohleos 
usually Involves coniplex mathematics and numerical techniques, which are 

nearly Iniposslble to treat without the aid of automated procedures. On this 
basis, a digital conqniter program based upon the Love-Reissner first order 

shell theory has been developed. Isotropic and kinematic hardening lavs 
are available using an orthptropic yield surface and including the 
Bauschinger effect. This program can perform a nonlinear gecmetric and 
material analysis of orthotropic thin shells of revolution, subjected 
to nonproport ional cyclic axlsymmetric distributed loaidlng or concentrated 
line loads, as well as thermal stralhs (Reference 1). Furthermore, a 
shell with artitrary boundary conditions, under loads which vary arbritrarily 
with position and under a temperature variation through the thickness, is 
tractable with this program. The shell can consist of any combination of 
the following geometric shapes: 

1) Ellipsoidal - spherical (offset from the axis of revolution allowed) 

2 ) Ogival - toroidal 

3 ) Modified ellipse shape 

4) Conical - circular plate 

5 ) Cylindrical 

6) General point input geometry 

7 ) Dummy geometry slot to be filled by the user 

8) Discrete ring 

The shell wall crossection can be a sheet, sandwich, or reinforced sheet or 
sandwich. The reinforcement can consist of rings and/or stringers, a waffle 
construction rotated at any angle to the principal coordinates, or an Isogrld 
cons t37uct ion. General stlffhess input options are also available. The rein- 
forcement material properties can differ from those of the main shell, and a 
temperature variation can cause different properties in the two face sheets of 
a sandwich shell. 
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The basic approach to the problem (Reference l) is to cut the structure into 
several shell regions. These regions need to be singly-connected shells^ 
and can only have line loads applied at their end points. There are no 
restrictions on geometry^ or uniform or thermal loads. The regions are 
further subdivided into several shell segments, each being free to have its 
own geometric shape, provided that the shape falls into one of the categories 
mentioned above. 

Stiffness matrices obtained for each segment, are coupled by standard matrix 
methods to obtain region stiffnesses, which, after being reduced in size, are 
in turn coupled to fom the total shell structure under analysis. Currently, 
the con^niter programs are sized to handle a structure coii^josed of up 
to 29 segments in each of 29 regions arbitrarily connected to each other. 

There is a limitation on the size of a shell segment, which is a consequence 
of the demand that boundary disturbances be felt throughout the segmen-^. This 
limitation is mathematically described in Section 2 (pages 2-35 to 2-37 ) as a 
length parameter. This parameter, however, is not reliable near the apex of 
any shell shape (c|) = O), and the segments needed in this region are actually 
much smaller than predicted by the parameter. A mathematical singularity 
occurs at the apex where r^ (the radius of revolution) becomes zero. It is 
this singularity which prevents the length parameter from being meaningful 
near the apex. Furthermore, the point (<(> = O) is not an acceptable input 
point of the program (except for the torus-ogive and offset ellipsoid), 
although any point outside a circle of infinitesimal radius is satisfactory. 

There is a considerable latitude in what can be done within each shell seg- 
ment. The thickness of any segment can be symmetrically tapered and it can 
contain up to 1^ points of discontinuity, provided that the segment center- 
line remains continuous and describable by a single shell geometry. A temp- 
erature distribution through the thickness can be specified at three points 
in a homogeneous shell, and h points in a shell of rigid core sandwich con- 
struction. The distribution is considered to be linear between these points. 
Thus, it is possible to approximate temperature distributions other than 
linear distributions. In the event of physically discontinuous shell center- 
‘lines, a kinematic link is available for use in the analysis- The link 
relates displacements across the discontinuity. This link may be used between 
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regions^ and between segments within a region. Discrete offset rings are 
also available for use within or between regions. 
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SECTION 2 


INPUT INFORMATION 


2.1 GENERAL NOTES 


The preceding section provides some insight into the capability of the pro- 
gram, and the potential that it might have for future use. If the program 
is applied Judiciously it can be an extremely powerful tool. The mech- 
anics of applying it should be clearly understood. With this in mind, the 
remaining section should be studied carefully. 

The required input data may be subdivided into three main parts, namely: 
geometric, topological (or coupling orientation) and Joint data (degree of 
freedom description for each Joint component). Each segment requires its 
own geometric configuration and numerical integration control. 

The output consists of stiffness coefficients for each shell segment and the 
actual symmetry of the coefficients is presented in a convenient form for 
a check on the accuracy of the integration through the segment. Region 
stiffnesses and their symmetry checks are also provided. Final stresses, 
displacements, and plastic strains are printed out for each shell segment 
at intervals along the segment as specified by the user of the program. 

The output will be further discussed in Section 3. 

The present program size is described in the table below. 

Table of Program Sizing 


I. 

Segments per region: 

29 

II. ■ 

Segment joints per region: 

30 

III. 

Regions: 

29 

IV. 

Region joints: 

30 
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Table of Program Sizing (continued) 


V. Number of points available 
per segment for specifying 
geometric or load data: 30 

VI, Number of points available 
through the thickness for 

specifying temperature data: ^ plus reinforcement ten^eratures. 

VII. Geometries: ellipsoid, sphere, offset ellipsoid, 

modified ellipsoid, ogive, toroid, 
cone, annular plate, cylinder, general 
geometry, ring, elastic support, dummy 
geometry. 

VIII. Wall cross-section options: single sheet, equal face sheet sand- 

wich, unequal face sheet sandwich, 
eccentric reinforcement (rings, 
stringers or both), waffle reinforce- 
ment rotated at an arbitrary angle to 
coordinate axes, isogrid reinforcement, 
arbitrary stiffness input. 


IX., 

Number of material property 



tables per submission: 

10 

X. 

Number of points per 



material property table: 

10 

XI. 

Hardening laws: 

kinematic, isotropic, perfect plasticity 

XII. 

Orthotropy options: 

isotropic or orthotropic sheet, iso- 


tropic or orthotropic sandwich, iso- 
tropic or orthotropic sandwich with 
different face sheet properties caused 
by themal gradients, isotropic or 
orthotropic sheet or sandwich reinforced 
by- different property- rings or different 
property stringers or both, isotropic or 
ortho tropic sheet or sandwich reinforced 
by a different property waffle system 
rotated by an angle p to coordinate axes, 
isotropic or orthotropic sheet or sand- 
wich reinforced by a different property 
general isogrid reinforcing system, 
arbitrary stiffness input options used to 
describe other configurations. 
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Table of Program Sizing (continued) 


XIII. Number of rings per region: 28 

XIV. Total number of rings at 

region joints in structure: 28 

Figure 2-1 shows the detailed option flow chart for the present programs. 


GENERAL NOTES - Idealizations 


Before discussing the specific card input order, it would be advantageous 
to introduce some general guidelines in the area of idealizations and 
topology. In many computer programs there is such an abundance of numerical 
computation, that minimizing numerical roundoff errors becomes as important 
as getting the final answers. In some cases the engineer can aid the 
program in this effort throught the use of judicious idealizations. Such 
a possibility exists in the STAES-2 programs, since many internal operations 
are involved with building and Inverting stiffness matrices. The object 
of the user therefore, should be to help the computer by avoiding the 
creation of ill-conditioned matrices at any step (see Reference 2). 
Physically, the way to achieve this end is to have all the segment stiff- 
ness matrices of the same order of magnitude. This will in turn produce 
region stiffness matrices which are of similar orders of magnitude, and 
min imize possible ill- conditioning in the total structure matrices. The 
user can help to achieve this end by sizing his segments in such a way 
so that no short stiff segment is contained alone in a region with all 
other long flexible segments, or that no region comprised of all short 
stiff segments exists in a structure whose other regions contain only 
long flexible segments. No accurate measure can be given on the relative 

stiffness or flexibility of segments allowed, and thus the best check is 
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Figure 2-1 Program Option Flow Chart 
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to if & BtniotuTfl 1 b ila cl^uilibriuAi luidBr tba applied loadlcigp In 
case of tba cboiiLLaear aaAlysic An internal equlUbriuiiL oorTBctloq 
ma^ be apetitled in the Input, synHBetry eheeta &f and region 

fltlffrteflfl mtrlcea aore uKful fw many reae^Dns^ biit %rill not noc^BBEarily 
alert a hebt to tll-condltioTLlngp 


In the use of regions ^ one other t/pa of accident must be avoided. This 
±B the creation of a sin^e rsgioo atr'icture vitb. both endB f±K.edj where - 
in no suitable bDut^dary condition ratrix can be formed, ThnSj, in the 
LiBe of reftlQct idealtaatioiia^ which are less ph/BlcaHy laeatiingfiil to a 
LLBer than para segnifint tdeallEatlonB , care Bhoald. be taken ao that all 
bKihdary con4itlons are not ^er«d out. To avoid thU probleiu.d to 
miniAixe progj^m fUnniEig timf, it i4 bMt to maxlDliX^ tho nUabOr Of 
region* in a stmcture^^ and Bininise the number of GogmootB per region* 
iniuB^ in small problems^ for best numerica.! efficiency^ there Bhonld. 
only be one segnient per regloii. 


In the solution of atatic axlaymmetric problaniB, torsional And non-^torBioiUal 
BtAten may bt unceupLe^d Thua if no toraionBl loadis ej^pt iji enfrnhi&atlon 
with tldn-tHSrsiOflftl load^p all torsional degr&M of freedom can he remoyedi 
When bath. typEB of loiuHA exist, th* nonlinear ajv&ly*is will be eoupled. 


Tin an inereoierLtal analy^lg oueh as the preBent^ the number of load incre- 
ment e is part of the user fipeeified "14eali^iAtiClb'^ eurroht anadyelfl 

uses each Lgad step m it one -step Uewton-Bji^pbaoi% cycle, Utus the cto^er of 
load Bteps necessary is diroctly proportioRAl to the tbonlinearity of tho 
material stress‘^^raij!i helLavlor^ lo genBral it may be necBBsary to run each 
problem twice t otELc* with large steps to d^t^rnine gross bel^avior mA then, 
after eoarparing the analysis dtrain ineranueaitB Uq tbA shell to the lAaterial 
property cnrYeB, with smaller stepE^ 



(jBfSSAl HDTE3 - flegtaj-bs 


SlciCfl solxirtloQ of a nodllEMflj' geotoetrlc find BiAteTlfiLl prableai raqiilrfifi 
th^ ccuBccut J_ve ejcecutloD of a EerieE of arialyEeg 3 tiie progrBon bo 
inadwrtwtly CUlt off due to ao Inauf^ljcleiit time EBtlnntej or the uaor 
my viflh tQ examine the Bolutiona of a i^ycliog or oanBecutlTa loa^i ng 
ptoblew befCJW executing further + In order to allow the user to utilize 
the gpoi dotfi. g^n^rated In awb ruus, a restart qptioQ has bean InoorporAtAd 
lnt& th* &!S 3 AB 3 - 2 p profiram* optloa is uaed as followa: 

After ttie Cfnnplate execution of any load step the ooUpleted l&ad 
Stop tlUPboT 1? printed, In a load cycling run both the crrerall 
ratep eounter (initially Etarted from zero) and the loeal step 
coon ter (initialized tn zero at every load change) *re jcrlEvtad- 
^ order to be ahle to restart^ ths tapes on units 11 5, l£j 13 and 
15 nniet he saved, (The first two Ma nfl^eSeary OElly if ttw run 
nontatcE dlEcrete rlnge, ) These are oiillled respertlvely EIHiCRll,, 
fiBPIASlJfc md 3HPLAai5^ 

S. Upon tcoelpt of M aborbed the oeceBEary restart ba^S c^an be 
Dbtaiaed froffl "the odPiJlfrted load ftep printout, If tha last 
cfverall load step conpleted wan oddj the neceEsary tapfifl are 
RDKIEEll and SKPIASlSi tti* other two saay ba rtloasod. If the 
last ever-all load fltep raapleted was even, the neceEEBry tapes are 
RES1SEB12 aad 3KPIA313; th* other two tiHiy bo rftltasedi (An 
exception to the above fcnmila ittvolvefl teBtartiag restArtfid runs: 
Tf the previous restairt started cm an even lead step tht above 
foimila is reverEsd. ] 

3. In ordar t& Initiate a restart rm the fbUcfiflns mat be done: 

a^ In a filngle load rua adjuM tbe Pro^am Ckhntrol Card as noted 
in Items 21 sfld SJ, 

b. Eft a load-time history or cycling run convert the load data 
by rtatvlng the cem^eted cyclaa from the deck» aud adjust 



the Program Control Card accordingly • Note the item about 
the Ramberg^Osgood parameter in the cycling input 
description on p 2-79* 

c« Operator instructions must be given to mount the SHPIAS 
restart tape on logical unit 13 9 and the RINGER restart 
tape (if necessaiy) on logical unit 11 (regardless of the 
numbers associated -with the restart tape names )» Scratch 
tapes will be mounted as usual on logical \inits 12 and 

d. The run will now proceed as noimal (even restarts of restarts)* 

4, Note: The restart tapes for any given run are not useable more 

than once. Upon the initiation of the computer run these tapes 
become scratch tapes and additional . infoicmiation (more load step 
data) will be written on ^hem. 

GENERAL NOTES - Data debugging 


The STARS programs have been provided with special separate data debugging 
packages called SATELLITE programs. In order to be able to debug as much 
of a given data deck as possible in one computer submission, the data is 
grouped by inserting special cards, termed "dash-separator cards", appro- 
priately. In order so that additional errors are not made by requiring 
insertion and removal of these cards, the STARS program has been coded to 
accept these dash-separator cards in the input. A dash- separator card is 
shown below : 





As can be seen, a minus • symbol is inserted straight across the computer card 
frcm column 1 through column 80 . 

Since the dash-separator cards subdivide the data deck, there exists the pos- 
sibility that a separated data block may be completely omitted (for example 
no kinanatic links in a structure). In this case one dash-separator card is 
also omitted . Under no circumstances can there exist, tvo adjacent dash- 
separator cards in a data deck . The SATELLITE programs are described in 
Reference 3. 



2*2 CAKD FORMTS 


The following different card formats are presently required by the STARS *-2P program. A full 
description and explanation of the information to be entered on the cards is presented in Section 2*3® 


Title Card 



(l) Alphameric Title 


16A4 


Program Control, Program Cycling Cards 


H 

Hi 

1 

II 

IB 

IBiBiB 

a$E 

a] 24 1 

m 

IBHIB 

|3S 

ilL 

~T3 

36 

7 3S 
38 

[41 1 ’43’ |45i l47| 

_4Pi 42. '44-, i46i 48 

49 1 SI 
so 

|53i 

1 S4 

H 

57 

58 

>9 lei 
60' 

63| |6S| 1671 68 7tl 


l®l 

@1 

© 





n 

M 

@ 


nj:-, 




IK 16^1 Iw| IQO fU 1/4 

/«! 1/6. jW 


(1) Rmber of regions 

(2) Total number of segments 

(3) Number of material property tables 

(4) Intermediate print clue 

(5) Number of load cycles 

(6) Number of load steps 

(7) Load step print increment 

(8) Load unit clue 

(9) Nonproportional load clue 

(10) New material table clue 

(11) Completed load steps in last cycle for restart 

(12) Total completed load step for restart 

(13) Blank 

(14) Rotation clue 

( 15 ) Rot at ionaJ. speed 


12 

13 

12 

12 

12 

F6.0 

F4.0 

12 

12 

12 

F6.0 

F6.0 

.12 

E14.7 




(1) Problem Identification 

(2) Equilibrium correction 

( 3 ) Blank 

( 4 ) Graphics clues 

(5) Graphics Print Cycle 


Material Property Table ID Card 


® 1 

(D 

l@l 




(1) 

Material Name 



(2) 

Blank 



(3) 

Table IJrpe 




Materials PrppOTty Table,. Ge(pmetry, Position, Crossection, Loading Cards 


\ 0 \ d) r 6) 

(l) Input table item (as many cards and fields are necessary) 

Region Introductory Card 


(1) Number of segments 

(2) Number of kinematic links 

(3) Number of rings 

(4) Title 




ST“S 


Topology Cards, (Region Joint Control Card) 




!#j| 

» !(■ 

09 




jsL 

£ 


F* 

M 

I 21 I j23| |?S! |Z7 
»' ' 22 ! 124; ' 2 (| a 

291 |3i: '33 |35| 

' 30 ! '32 34 ; !3$ 

37 

|39; l4i! 43’ . 45 ’ ‘47 t49| 'si! [s3: Iss] !s7j » III 

39! <40, ‘42 44 . 41. 4» SO S7 ' S4. S4i SI 10 

]S3 1$) Il7 

12 M in M 

IIJ 

71 

A 

73 

74 

7si Inf 179 j 

\n: \n\ 

a> 

0 1 

d) 



(1) Region or segment number, (number Of joints) I5 

(2) Beginning joint, (number of rings) I5 

(3) End joint, (number of links) I5 


Segment Identification Card 


foil 

i 02 


y 


1 01 

0 

Inl 

l»> 12 

0 

l»| 

u| <11 

0 

Iw 

III 

l 2 li | 23 | | 2 S! |27 
20 ' ' 22 ' ! 24 | s' 

B 

”1 

”30 

IV <33 
32 


a 

0 

| 39 i ! 4 li ’ 43 ’ i 4 S: ‘ 47 V i 49 | < 5 l! l 5 s|~[wr 

31 ! I 40 . '42 44 ‘ 41 . 41 ' SO S 2 ' $ 4 . S|! '.SI 

0 

m. 

“I« 

e 

Q 

d ^!»2 

!«! 

1 « 


73 | | 75 < 1771 | 79 i 
I 24 I 174 , ;« 

\® 

m 










d 

) 






(1) Segment code F2.0 A1 

(2) Number of plastic layers 12 

(3) Title 16A4 


MAGIC Integration Card 


t i| i(B; io$[ 

i«i ,« l l« 


07 1 09 U 13 

| 0 I| |l 0 | 12 | }U 


IS 17 19i 12! 23 2S 

»{ |I«r |20{ 22 24 


© I @ 

(1) Print Interval 

(2) A«curacy control 

( 3 ) Integration interval 

(4) Blank 

(5) Step control 


3 t I • 133 I 3 S 


[43 ‘ 45 ; i47' 

j j40} ,42 ig ^4t; |4» 


I *3 
jS2 ; ;S4 


Mi {57 
ISO 


& 


«I~| 7 l| I73T fcl 

ml 72f 74 I 70 


E14.1 
El4. 1 
El4. 1 

F2.0 


Special Geometry Card for Point Input Shape 



® I © (D \ 0 


0 0 0 


(1) Glues 

(2) Blank 

( 3 ) Number of Table points 










f{T-Z 



(1) Joint Number 

(2) Boundary conditions 

(3) Axis rotation angle 


12 

F2.0 

E14,1 
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(1) N-umber of joint loads 

(2) Title 


Joint Load Cards 



(1) Load condition nmber I 5 

( 2 ) Row Identification I 5 

(3) Applied load EIU .7 





2.3 DETAIIED ORDER OF INPUT (See Figures 2-1, 2-2) 


GENERAL INTRODUCTORY CARDS 

Column 

Format 

1. 

Title Card 




A, 

Alphameric title (submission description) 

1-64 

16A4 

2. 

Program Control Card 




A. 

Number of regions to be coupled 
(Max. = 29) 

1-2 

12 


B. 

Total number of segments 

(Max. = 29 X 29 = 84l) 

3-5 

13 


C. 

Number of Material Property Tables 
(Max. = 10) 

6-7 

12 


D. 

Intermediate print clue 

8-9 

12 


If this clue is set to unity, stiffness 
matrices and synmetry checks -will be printed 
for all segments and regions at the specified 
load step print interval (item G below). To 
delete this intermediate print the clue should 
be set to zero. 

E. Number of load changes 10-11 12 

Nunober of load patterns or cycles to be applied 
in a cycling load ^alysis (see p 2-76 ). 

F. Humber of load steps 12-17 f6.0 

Number of loading increments to be made for 
the first loading cycle. 
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Figure 2-2 Data Sequence 
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Column 


Format 


G. Load step print increment 18-21 

The interval niimber of load increments at vhich 
printout is desired. (Ex: print every 20 steps). 

H. Load unit clue 22-^3 

This clue is zero (O) for a single loading 
analysis. If there is a change of load 
distribution with time (sequential loading) 
this clue is two (2). For further information 
and load cycling see the input for the Losid Cycling 
Program Control Card on p 2-?6 . 

I. Local cycle load step count for restart 28-33 

If the restart submission is concemed with load 
amplitude changes or cycling, this aatry is the 
local count of load steps completed with the last 
load ing amplitude input. If there is no load amp- 
litude changes item J below and the present item 
are equal for restart. (See p 2-7 for restart 
information. ) 

J. Total load step count for restart 3^-'!’39 

In the event that the submission is a restart 
problem, this entry provides the total count 
of previously conqpleted load steps. 

K. Shell rotation clue 57-58 

If the shell is rotating at an angular velocity 
about its centerline, the input clue is uni^. 

If there is no rotation the clue is zero . 

L. Angular Velocity in rad. /sec. 59-72 

The maximum an5)litude of angular velocity of 
the rotating shell. This velocity will be 
attained in a series of steps in a nonlinear 
analysis analogous to a mechanical or thermal 
loading. The step size has been defined in 
item F of this card. 


FU.O 


32 


F6.0 


F6.0 


32 


E14.7 
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Column 


Format 


3. Analysis Option Card 

A. Analysis Description Clue 1-4 A4 

There are three possible analysis formulations. 

If the problem is to be a smal 1 deflection 
plastic analysis, then the clue word is PIAS . 

If the problem is to be a large deflection 
elastic analysis then the clue word is MLIN . 

If the problem is to be a large deflection 
inelastic analysis then the clue word is KLPL . 

B. Equilibrium Correction Clue 5-10 l6 

Equilibrium correction loads can be used with 

a.l 1 three analysis options above. To utilize 
this option set the clue to unity . To eliminate 
equilibrium corrections set the clue to zero . 

It is recommended that equilibrium corrections 
always be included. 

C. Graphics Control Clues 

The following nine (9) data fields determine the 
plotting required for the run. 


a. Circumferential displacement clue 21 II 

Clue = 1 plot u 

Clue = 0 or blank, do not plot u 

b. Meridional displacement clue 22 II 

Clue =1, plot V 

Clue = 0 or blank, do not plot v 

c. Noimial displacement clue 23 II 

Clue = 0, plot w 

Clue = 0 or blank, do not plot w 

d. Hoop stress (inner face) clue 24 II 

Clue = 1, plot CTq 

Clue = 0 or blank, do not plot CTq . 

' o in 

e. Meridional stress (inner face) clue 25 II 

Clue = 1, plot CT . 

’ 9 in 

Clue = or blank, do not plot a 

’ cp in 
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Format 


f. In-plane shear stress (inner face) clue 
Clue =■ 1, plot 

Clue = 0 or blank, do not plot t „ . 

cp9 in 

g. Hoop stress (outer face) clue 
Clue =1, plot ci 

0 out 

Clue = 0 or blank, do not plot Oq 

h. Meridional stress (outer face) clue 
Clue = 1, plot CT 

Cp OlXv 

Clue = 0 or blank, do not plot a 

cp out 

i. In-plane shear stress (outer face) clue 

Clue = 1, plot T „ . 

’ cp0 out 

Clue = 0 or blank, do not plot t „ . 

cp6 out 


D. Graphics Print Cycle 

Graphical output will be provided in accordance 
with the Graphics Control Clues above at every 
’X’ load step as input here. 

miERIAL PROPERTY TABLES (Max. = 10 sets) 


As many sets of these cards are used (^0) as there are 
different material property segments in the structure 
to be analyzed. These tables will 'be used to obtain 
the thermal variation of material properties if thermal 
loadings exist. Thus the range of temperature in this 
table should be greater than that of the thermal loads. 
If no thermal loads exist, the values given in the first 
column of this table will be used, and the rest of the 
table can be left blank. If there are thermal loads, 
the range of the table is to be considered as that 
between the second and tenth columns. 

Since apex segments should not be allowed to go plastic, 
separate tables with fictitiously high yield stresses 
should be prepared for this case (see p. 2-35 ). 


26 


27 


28 


29 


32-35 


II 


n 


II 


II 


fU.o 
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MATERIAL PROPERTY TABLES (continued) 


ColtJinn 


Format 


1. Identification Card 

A. Material Title (Alphameric) 1-^ 

Any name can be made up as long as it is con- 
sistently used on the segment cards to which 
it refers. The same name cannot appear on 
more than one (l) table. 

B. Type of Table 11-1^ 

One of several possible alphameric clues is 
written here. These clues serve to size the 
number of cards in the property table, and 
define which properties belong on which card. 

The possible clues are: 

ISOT 

ORTH 

STIF 

Their definitions are provided in Item 2 below. 

2. Material Property Cards 

The material property cards below are given depending 
upon which table type clue is used. If the table 
type clue is "iSOT" (isotropic table): 

A. Temperature values (5 values per card; 2 cards) 

These are the temperatures at which the values 
of material properties will be given. The first 
value in the table must always be the room or 
stress-free temperature, since the material 
properties in only the first column of the table 
will be used in an analysis involving no thermal 
load-. The values of temperature 'in table col- 
umns 2 through 10 must be in algebraically 
increasing order. 

B. Values of Young’s Modulus at the given temperatures. 

(5 values per card; 2 cards) 

C. Values of Poisson’s Ratio at the given temperatures. 

(5 values per card; 2 cards) 

D. Values of the thermal coefficient of expansion 
at the given temperatures. 

(5 values per'^card; 2 cards) 


Ah 

Ah 


5E14.7 


5E1^.7 

5E1^.7 

5E1^.7 
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E. Values of the Ramb erg-0 sgood ”s" (see Ref. l) 
at the given ten^eratures. (5 values per card; 

2 cards). 

F. Values of the Ramberg-Osgood "h” (see Ref. 1) 
at the given ten^jeratures. (5 values per card; 

2 cards). 

G. Values of the yield stress at the given ten^iera- 
tures. (5 values per card; 2 cards) 

Mote ; If the hardening law for the segment is 
perfect plasticity (PERF on the Master Clue Card) 
the items E and F. above are zero . 

If the table type clue is "ORTH" (orthotropic table); 

A. Temperature values 

(5 values per card; 2 cards) 

These are the temperatures at which the values of 
material properties will be given. 

B. Values of Young's Modulus in the 9 direction (Eg) 
at the given temperatures. 

(5 values per card; 2 cards) 

C. Values of Young's Modulus in the <|> direction (E ) 
at the given temperatures. 

(5 values per card; 2 cards) 

D. Values of the Poisson's Ratio v at the given 

temperatures. ^ 

Orthotropic identity defined as = ^cp'cpS' 

(5 values per card; 2 cards) 

E. V^ues of the thermal coefficient of expansion in 
the 9 direction (oCq) at the given temperatures. 

(5 values per card; 2 cards) 

F. Values of the thermal coefficient of expansion in 
the <() direction («^) at the given temperatures. 

(5 values per card; 2 cards) 

G. Values of the Shear Modulus G at the given 
temperatures. 

(5 values per card; 2 cards) 


F_ormat 

5E14.7 

5E14.7 

5E14.7 


^Elh.7 

5E14.7 

5E1^.7 

5Eli^.7 

5E14.7 

5E14.7 

5Eli»-.7 
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Values of the Ramberg - Osgood "s’* in the cp direction 
(see Ref* l) at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

I* Values of the Ramberg - Osgood "n” in the cp direction 
(see Ref* l) at the given tenperatures* 

(5 values per card; 2 cards) 

5E14,7 

J* Values of the yield stress in the cp direction at the 
given temperatures* 

(5 values per card; 2 cards) 

5E1U.7 

K* Values of the yield stress in the 9 direction at the 
given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

L* Values of the yield stress in the normal ( q ) direction 
at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

M* Values of the shear yield stress in the cp-G plane at the 
given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

IT* Values of the Ramberg - Osgood "s" in the© direction 
(see Ref* l) at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

0* Values of the Ramberg - Osgood ^!n" in the e direction 
(see Ref* l) at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

P* Values of the Ramberg « Osgood "s" in the cpG directions 
(see Ref* l) at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

Q. Values of the Ramberg - Osgood *!n” in the cp0 directions 
(see Ref, l) at the given temperatures* 

(5 values per card; 2 cards) 

5E14.7 

Rote; If the hardening law for the segment is perfect 
plasticity (PEKF on the Master Clue Card) items H, I, 
and N through Q are zero. 



2-^23 



Colxamn 


If the table type is ”STIF” (table to be used for 
reinforced shelJLs) ; 

A.-Q. The values in these locations are the same as 

those above for the "ORTH" clue case, and refer 
to the basic shell. 

R. Values of ring Young ^s Modulus (E^^) at the 
given temperatures. 

(5 values per card; 2 cards) 

S* Values of stringer Young’s Modulus (E ) at the 
given temperatures. 

(5 values per card; 2 cards) 

T. . Values of ring thermal coefficient of expansion 

(o^) at the given temperatures. 

(5 values per card; 2 cards) 

U. Values of stringer thermal coefficient of 

expansion (a ) at the given temperatures, 
s 

(5 values per card; 2 cards) 

V. Values of the Ramberg - Osgood "s " in the 
rings (see Ref, l) at the given ten 5 )eratures . 

(5 values per card; 2 cards) 

W. Values of the Ramberg - Osgood "h" in the 
rings (see Ref. l) at the given ten 5 )erat\u‘es . 

(5 values per card; 2 cards) 

X. Values of the ring yield stress at the given 
ten 5 )eratures . 

(5 values per card; 2 cards) 

Y. Values of the Ramberg - Osgood "s” in the stringers 
(see Ref. l) at the given temperat\ires . 

(5 values per card; 2 cards) 

Z. Values of the Ramberg - Osgood "n" in the stringers 
(see Ref. l) at the given temperatures. 

(5 values per card; 2 cards) 

a» Values of the stringer yield stress at the given 
temperatures . 

(5 values per card; 2 cards) 

Nfotes ; If the hardening law for the segment is 
perfect plasticity (PERF on the Master Clue Card) 
items H, I, N through Q, V, W, Y, and Z are zero. 


Format 

5E1^^.7 

5E1^^.7 

5E14.7 

5E1^.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 
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In a rotated waffle or isogrid constructions, items 
R and S, T, and U, V and Y, W and Z, X and ry refer 
to the grid directions and are respectively identical. 

D-A-S-H S-E-P-A-R-A-T-O-E C-A-R-D minus in l-80 

(See General Notes - Data Debugging) 


REGION INTRODUCTORY CARDS 


These two cards are placed at the beginning of each region 
data information. Each region contains the following data 
set (see Figure 2-2): a) Two region introductory cards; 

bV data cards for each segment within the region; 
c) ring cards describing the discrete rings within the 
region, if any; and d) kinematic link cards describing 
the kinematic links within the region, if any. 

1. Identification Card 

A. Number of segments within the region (< 29 ) 1-2 

B. Number of kinematic links between segments 3-^ 

within the region. 

C. Number of discrete rings between segments 5-6 

within the region. 

D. Any alphameric information (region description) 7-70 

2. Topology Card (Coupling Orientation) 

A. Region Number 1-5 

Number of the region under consideration. 

B. Joint (i) 6-10 

Joint associated with i (beginning) end of 

the region. 

C. Joint (J) 11-15 

Joint associated with J (ending) end of 

the region. 

There is no coordinate flow in regions (unless 
1 region = 1 segment), such as that shown for 
segments in Figures 2-3 to 2-9<> However, the 
start joint of a region must match with 1 in 
segment numbering, and the end Joint must 
match with the highest segment Joint number 
in the region (see Figure 2-12 and page 2-57). 

If 1 region = 1 segment^ the segment topology 
car^ will be a dummy card [ 1 12 ] (see page 
2-57). 


12 

12 

12 

i6a4 

15 

15 

15 
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Figure 2 - 3 - Typical Shell Segment 
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Figure 2-4- a. Ellipsoid 











1 



I 


n = -1/2 b/a = 0.639 
n-^-1 b/a -►0.618 


Figure 2-5. Modified Ellipsoid 
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z 
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i in j Specify (f> = 0 


J 

Inverted Cone 
Specify d) < 0 



Figure 2-8. Cone 






+z 



Specify Tq 


Figure 2 - 9 * Cylinder 
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SEGMENT CAEDS 


Column 


Format 


This sequence of cards is repeated for each segment 
within the region. 

1, Identification Card 

A. Segment identification code (see Figures 2-h to 2-9) 


1. 

Ellipsoidal or spherical shell 

Code 

= 11 

1-2 

K.O 

2. 

Modified ellipse shape 

Code 

= 12 

1-2 

F2.0 

3 . 

Ogival - Toroidal 

Code 

= 13 

1-2 

F2.0 

4. 

General Geometry 
(see Figure 2-7)* 

Code 

= 14a or 
i4b 

1-3 

F2.0A1 

5 . 

Dummy geometry slot 
(<{i coordinate) 

Code 

= 15 

1-2 

F2.0 

6. 

Conical - Circular Plate 

Code 

= 21 

1-2 

F2.0 

The plate is treated as a cone with 

zero 

angle. 



7 . 

Cylindrical shell 

Code 

= 31 

1-2 

F2.0 

Number of layers through the thickness 


4-5 

12 


This is the number of layers through the thickness 
for plastic strain integration. In a single sheet 
construction the maximum number of layers is 20. 

In a sandwich construction both face sheets are 
divided into the same number of layers and this is 
the required input. The maximum ntunber of layers 
for the sandwich face sheet is 8. 

C. Any alphameric information (segment description) 6-67 16 AU 

2. "magic" Control Card 

A. , Intesrval at which final answers are to be printed l-l4 El^.l 

(in radians or inches). 

The (^-coordinate is defined for all geometric 
shapes except the cylinder, cone and plate, 
for which the s coordinate is used. Figures 
2-k through 2-9 describe these coordinates 
for each shape. 

B. Difference 15-28 El4.1 

The value recommended depends upon the 
computer used. For eight figure accuracy 
computers it is 1.0 E-6; for the IBM 36 O 
it is 1.0 E-^. 

C. Integration interval 29-h2 El4..* 

The Runge Kutta numerical integration procedure 
is ■ substa.nt.-i ally more accurate than finite dif- 
ferencing. An interval of (.01 to.03) x segment 
size (in radians or inches) is recommended in eigen- 
value analysis so as to be able to represent eigen- 
vectors with high wrinkling in the segment interior. 

(in using a 30 point segment -bable (see p.2-4l)there 
should be at least 30 integration steps.) 
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D, Apex Clue 46-49 A4 

In the STARS program an apex segment is liable to 
have fictitiously high stresses due to the sing- 
ularity effect of r^ -» 0, Thus such a segment 

should not be allowed to go plastic ^ since it may 
prematurely stop the rim. Therefore a separate 
material property table should be prepared containing 
fictitiously high yield stresses. The word APEX 
should also be entered here as a reminder ^ and will 
trigger a reminder message in the output. 

E. Delta 71-72 F2.0 

For a fixed-step integration^ Delta = 0. 


This card controls the Runge-Kutta numerical 
integration scheme. The suggested values above 
yield accurate results for a fixed-step integration 
method. 


Calculation of Segment Length 


There is a restriction on the length of the shell segments. Physically^ 
the restriction demands that boundary disturbances at one edge be 
distinctly felt at the other edge. This is a consequence of using a 
numerical integration procedure. Since the segment stiffness matrices 
must be symmetric, the calculations involved in obtaining each matrix element 


must be such that a computer round off error never becomes prominent. Limit- 
ing the segment length insures satisfaction of this criterion. This length 
is a function of both geometric shape and segment location within a specific 
geometry. One of the limiting factors is that the ratio of the radii of 


revolution at the initial and final points of a segment be greater than one 


hundredth and less than one hundred . Thus < 


<100 where: 



This requires smaller segments than will normally be predicted by formula in 
the area of an apex. In addition, note that (cf) = O) is not an acceptable input 
point (except for the torus -ogive or offset ellipsoid). Also the apex segment 
should not be allowed to go plastic by adjusting the corresponding material 
property table . 
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For a cylinder^ the segment length parameter^ 


A=(1-i*y)^6As 

should be held to about k.O. In this expression^ ”7” is a non-linear param- 
eter. For homogeneous shells i 



It is zero for a linear problem. 

"P” is a measure of the rate of .decay of a disturbance in the shell. 
”As” is the meridional length. 


r 

o 


sin (j) 



Approximate formulas can be obtained for near cylindrical regions of generally 
curved surfaces. The length parameter , 



(1 + y )^ 


X A<p 


should be held to about U.O. In this expression ”y” has the same definition 
as in the cylinder case. 
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”X” is a measure of the rate of decay of a disturbance in the shell* 


^ is the angle intercepted by a meridional arc length ”ds"* 
The values of and As for various shell geometries are given below; 



-3 

The minimum allowable segment length is 1 x 10 
(inches or radians according to segment sizing). 


3. Geometric Description Card 


Ellipsoid and sphere (Figure 2-^) 

Column 

Format 

1. Semi-axis perpendicular to Z-direction (a) 

l-l4 

E14.1 

2. Ratio of semi-axis in the Z-direction (h) 
to (a), ^ - (1) 

15-28 

E14.1 

3. C = offset distance (±) (C = 0.0 if no offset) 
Modified ellipse shape (Figure 2-5) 

29-42 

E14.1 

1. -Axis ratio coefficient (n) 

l-l4 

E14.1 

2. Semi-axis perpendicular to Z-direction (a) 
Ogive (Figure 2-6) 

15-28 

E14.1 

1. = radius 

l-l4 

E14.1 

2. C = offset distance (±) 
General Geometry (Figure 2-7) 

15-28 

E14.1 

1. Number of pairs of Z versus r^ points (<l4) 

1-2 

12 

2. Z versus r^ points, in pairs, starting 

3-72 

FIO.O 

with Z value (7 values per card, including 
first card, for up to 4 cards total). 

1-70 

FIO.O 


In the input table the first Z value is 
taken as =0, and furthermore Z. > 0 

(i = 2 - l4) (see Figure 2-7). 
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E* 


Cone £-0) 

1- Ane3.e + in Todians flat plate, + = 0). 

Eoep la mind that thin + is a constant f&r 
a given, wn* ir'd ahaold not be confusad 
^Ith ttlfl! 4" on eard set 5- 


1-lb 


F4 Cylinder (Figure 2-9) 

li EadlLia 1-14 

4. Master Clue Card 

TMg Card docitaius a Eerie s &f -eluoa which da- 
te mtiirw? the program a^d tahlo- locations to be 
used for the aegmetit defloribed, Fbr a 

raaater flow ehart of oliiea and options in the 
progr^ffn set Figure 2-1. 

Ar Material Table Typo Clue 1-4 

This clue derives the type of material property 
table to be CKpected for tjie assent- Thts^ aS 
well clS the following cluE detenniues the mate- 
rial propertlea th^t will be used in the strud-- 
tural anBlycis for the aegment, Thna theae two 
oli^ea ahould match the ti;ro «lues used on tho 
Idetitlflnatioti card of the corre ponding mate- 
rial property table. As mentioned before oti 
pageS--21 , the three posSihilitiefl are: 

ISOT 

OETK 

STIT 

B- Maitcrlal Title 12-14 

This name* Ehould be the same aa the name -which 
appear a on the material property table which 
containa the properti-ee td be utilla&cd for 
thia aegment. 

C. Sheet Clue 21-24 

This cine InrorfDS tho pro^am be to What kind 
of ahfrll Wall eroBBflc-tion to e^act. If the 
shell is of a ingle aheet c&nstruatiqn^ th* clue 
to be used isi SITfSi If the shell wall la an 
actual- aiae face aheet e{Wi4c^oh^ the clue to be 
uaed te: EQUA+ If the ahell wall ia a aand- 

wleh but the faoe aheeta are not equals tho 
clue to be use^ is; FinallyV if the 

shell segment la reitif^arded by rings, stringers j. a 
waffle j or Eih leogrid^ the cLiie to be used 
is: BIiAH . 


2*36 




EL4.1 


A4 


At 


A4 





Fcirwflt 


D. Reinfoi^eraGnt Clue 3l*3^ ^ 

This cJ_Lie deeoidles the tjT^ rtLnroffcePbiHt 
tbit ifi pPcaent 6fi the- ebcH- tf tbe Ebell is 
purely of single sheet or o^usl or 
face Eheet honey oolat coKifitTuetloci (no Tetnforclo^) , 
tht tlufr to be a&ed. Is: THIC. If tba re-triforce- 

nent eoniileta of rir\$& or ptrltigara or hothj 
located along the coordinate axes (S afld ^ or a)ji 
three a^t poa^iblB defpenfllng upon the 

baeic shell wall conetriAottoni If the h*slc 

1e b single sheets the clue to 
he iJE-ed is STU x If the basic shell vail i$ an 
equal- else fane ah«t sandwich j tbo clUe to he 

£□ ETLB d If the basic shell VbII is a Ssnd- 
’jlch but the fa-na shaeta bt^ ntiequalj tb^ to 
ho uied is ET13 i If the rEinfOrceaent oonslsta 
of a vaffle uHlch la rotated at ah arhltraj^y angle 
D from tbe meTidiondl -erlfl the following three 
eluea art pofiEiblc depending uppn the basic shell 
wall construction* If the haste wall e&natruction 
is a Single sheets the due to be used is HWAl-^ 

If the baaic ehell wall Is equal-Si-^t fftpe 
Sheet BaiidWidh^ the clua to be used is gtfAg. If 
the basic shall wall is a scjidwlch hut bha face 
shoets unequal j the clue to bo used £a FtfAj . 

If the rainfoToerjent ootielate ox\ isogrid eon- 
EtmetioR of gef>eral angle from tha meridLiOn al 
■AXia fnervtally ^ = 30'") the fallowing three clue? 
ore poasible depending upon the basic eheli wall 
most ruction- If the basic wall Cou£tr«nitiDn is 
4 single sh^e-tj the cHw to be used is I SGI , If 
the basic shell Wall is an equal-size face sheet 
sandwich j tbe clue to he used is If the 

basic shell vail is a Sandwich hwt the face sheets 
are unequal j the cluo to be used is lEjH - Two 
other clues are av^tLablej naffnsly S-TIO and 
if the user wlEhee to input hid own stiffnese 
constants- OTnStajltS may represent BTiy 

waU, construction -as long an the baBlo Hdo|I^''s 
L aws used with the Clnta are appropriate to 
describe the construction to be oonEidered- The 
Koote's Laws used with thepe cluas are gtv^n 
under the- dfrftcriptlon of segment card set 6 . if* 
plitJtie M&lyB±E is ellawed with the hast two 
Clu*B einOe a geametrin descrlpticfl of thfl croSdettica 
la not input. 
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Fottl 


Note ; The reinforcement described in the seg- 
ment cards is closely spaced reinforcement which 
will be smeared over the segment • Discrete rings 
at segment ends in a region are described at the 
end of all segment data for that region, 

E. Thermal Clue 4l-44 A4 

This clue describes the type of thermal problem 
which exists in the segment. The user is reminded 
that in a cycling run this clue is constant for all 
cycles 3 requiring zero thermal loads to be input 
-if necessary. If there is no thermal load on the 
segment, the clue to be used is NOTH , If the 
thermal loading on the segment is of general, 
standard type, that is if there is variation 
of temperature through the thickness as well as 
in the coordinate directions, the clue to be 
used is THST , If the thermal load is such that 
the variation is all in the coordinate directions, 
and there is no thermal variation through the 
thickness, the clue to be used is THCN , The last 
clue concerns a shell which is inhomogeneous in 
the meridional direction. This is not really a 
thermal problem at all, but merely a manipulation 
of the material property tables. If a structure 
has a wide variation in material properties in the 
meridional direction, without this last option one 
must take short segments of constant properties 
for analysis. With this option, however, the 
property variation is placed in the material 
property table, and expressed on the segment as. 
a function of temperature. No thermal loads are 
calculated, however, and the temperatures are 
only used to interpolate for material properties 
as integration is progressing along the segment. 

Thus continual variation of properties in the 
meridional direction is accommodated. The clue 
for this option is THIN. 

F» Hardening Law Clue 51-54 A4 

This clue determines the hardening law to be used for 
the segment. If isotropic hardening is to be used 
the clue is ISOT . If kinematic hardening is to be 
used the clue is KINE . If perfect plasticity is 
to be used the clue is PERF . 

Note ; ISOT hardening law is not allowed with ORTH 
material (item A on this cardJT^ 
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Column 


Format 


G • Table control - Number of points in each of the 71-72 12 

following tables* 

This can vary from 2 to 30 depending upon the 
shell geometry and loading. For a linearly 
varying geometry and/or loading only 2 input 
points would be required. These two points 
would be the end points. For more general 
loading and/or geometry a larger number of 
points are required. Ih particular, each 
abrupt change is specified by two points. 

One should use as many points as necessary 
(up to 30) in order to completely describe 
the problem, rather than using very short 
segments. 

I;. Table of <\> or s Values 

A. Initial, intermediate and final values of 5E1^.7 

<() or s. Each point requires l4 columns on a 
card and thus there can -be 5 values per card 
and up to 6 cards to make a total of up to 
30 points. The points to be specified are 
the beginning point of the segment, any point 
of discontinuity, and the end point of the 
segment. The input must be consistent with 
item G of the previous card. 

6. Table of Wall Crossection Geometry 

The contents of these cards (up to 6 cards per 
item below) are dependent upon the clues reg- 
istered on the Master Clue Card. If the shell 
to be described contains no reinforcing, the 
pertinent clue is item ^C, the Sheet Clue. The 
geometry is input and the stiffnesses are cal- 
culated by the program (see Figure 2-10). The 
input is presented below as a function of the 
Sheet Clue. 

If the Sheet Clue is SING (single sheet con- 
struction) : 

A. Initial, intemediate and final values of wall 5El4.7 

thickness (h. ) at points defined by table of 
or s values. 


\ 
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Fomat 


If the Sheet Clue is BQUA (equal-size face sheet 
sandwich) : 


A. Initial, intermediate and final values of face 
sheet thickness (h. ) at points defined by table 
of ♦ or s values. ^ 

5Eli^.7 

B. Initial, intermediate and final values of core 
thickness (t) at points defined by table of <|) 
or s values. 

5E14.7 

If the Sheet Clue is UWEQ (unequal-size face sheet 
sandwich) : 


A. Initial, intermediate and final values of inner 
face sheet thickness (h. ) at points defined by 
table of (|> or s values.^ 

5E14.7 

B. Initial, intermediate and final values of core 
thickness (t) at points defined by table of <|) 
or s values. 

5E1U.7 

C. Initial, intermediate and final values of outer 
face sheet thickness (h ) at points defined by 
table of (|) or s values. 

5E14.7 

If the shell is reinforced, the Sheet Clue will be 
BLAN. In this case it is the following, or Reinforce- 
ment Clue (item 4 d) which will determine the contents 
of card series 6 . For the reinforcement cases the 
geometry can be complex and varied, since all types 
of reinforcing are to be included. The reinforced 
shell input is presented below as a function of the 
Reinforcement Clue. 


If the Reinforcement Clue is STll (single sheet 
reinforced by rings and/or stringers): 


A. Initial, intermediate and final values of the 
torsional stiffness in the <|i direction (GJ ) at 
points defined by table of <|) or s values. 

5E14.7 

B. Initial, inteimediate and final values of the 
torsional stiffness in the 6 direction (GJ«) at 
points defined by table of (j> or s values. 

5Elij-.7 

C. Initial, intermediate and final values of 

stringer area (A ) at points defined by table of 
41 or s values. 

5E14.7 

D. Initial, intermediate and final values of ring 
area (Ag) at points defined by table of 4 or s 
values. 

^Elk.7 

E. Initial, intermediate and final values of 

stringer eccentricity (measured inwards from 
base shell centroid as positive) at points 
defined by table of 4 or s values. 

5E14.7 
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F. Initial^ intermediate and final values of ring 
eccentricity (measured inwards from base shell 
centroid as positive) at points defined by 
table of (}> or s values* 

G. Initial^ intermediate and final values of 
stringer moment of inertia (about base shell 
centroidal axis) at points defined by table 
of oi* s values. 

H. Initial, intermediate and final values of ring 
moment of inertia (about base shell centroidal 
axis) at points defined by table of or s 
values. 

I. Initial, intermediate and final values of 
stringer spacing at points defined by table 

of (|> or s values. (Do not set to zero if no stringers.) 

J. Initial, intermediate and final values of ring 
spacing at points defined by table of 4) or s 
values. (Do not set to zero, if no rings.) 

K. Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of 4> or s values. 

If the Reinforcement Clue is ST12 (equal face sheet 

sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the STll clue above. 

K. Initial, intermediate and final values of base 
shell face sheet thickness (h. ) at points de- 
fined by table of or s values. 

L. Initial, intermediate and final values of base 
shell core thickness (t) at points defined by 
table of 4> or s values. • 

If the Reinforcement Clue is ST13 (unequal face 

sheet sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the STll clue above. 

K. Initial, intermediate and final values of base 
shell inner face sheet thickness (h.) at points 
defined by table of 4) or s values. 

L. Initial, intermediate and final values of base 
shell core thickness (t) at points defined by 
table of 4> or s values - 


Fomat 

5Elij-.7 

5KL^+.7 

5E14.7 

5E14,7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 
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M, Initial^ intermediate and final values of base 5El4,7 

shell outer face sheet thickness (h ) at points 
defined by table of <() or s values. ^ 

If the Reinforcement Clue is RWAl (single sheet 
reinforced by a waffle rotated at an arbitrary angle 
from the meridional direction): 

A. Initial, intermediate and final values of waffle 5E1^.7 

grid area at points defined by table of 4) or s 

values. 

B. Initial, intermediate and final values of waffle 5E1^.7 

grid eccentricity (measured inwards from base 

shell centroid as positive) at points defined 
by table of (j) or s values. 

C. Initial, intermediate and final values of waffle ;>El4.7 

grid moment of inertia (about base shell centroidal 

axis) at points defined by table of 4> or s values. 

D. Initial, intermediate and final values of waffle 5E14.7 

grid spacing at points defined by table of or 
s values. 

E. Initial, intermediate and final values of waffle 
grid rotation angle, P, (in radians from the 
meridional direction) at points defined by table 
of or s values. 

F. Initial, intermediate and final values of extreme 
distance from base shell centroid to waffle outer 
edge (± value, positive inwards), at points 
defined by table of <p or s values. 

G* Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of {() or s values. 

If the Reinforcing Clue is RWA2 (equal face sheet 
sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction): 

A. through F. The items contained on these cards are 
those described for the RWAl clue above. 

G- and H. The items contained on these cards are 2 sets of 5El4.7 

those described for the ST12 clue above as items 
K. and L. 

If the Reinforcing Clue is RWA3 (unequal face sheet 
sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction): 

A. through F. The items contained on these cards are g sets of 5El4'.7 

those described for the RWAl clue above. 

6. through I. The items contained on these cards are 3 sets of 5E1^.7 

those described for the ST13 clue above as items 
K. through M. 


5E1^.7 


5E14.7 


5E14.7 


6 sets of 5E1^.7 
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If the Reinforcing Clue is ISGl (single sheet 
reinforced hv a general angle isogrid construction); 

A. through D. The items contained on these cards 
are identical to those described for the RWAl 
clue above ^ but with reference to the isogrid. 

E. Initial, intermediate and final values of the 
isogrid angle, p, (in radians from the 
meridional direction, see Reference ^ 

Appendix A). For the normal isogrid, the 
angle is 30° (input in radians). 

F. Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of <() or s values.^ 

If the Reinforcing Clue is ISG2 (equal face sheet 
sandwich reinforced by a general angle isogrid 
construction) ; 

A. through E. The items contained on these cards 
are those described for the ISGl clue above. 

F. and G. The items contained on these cards are 
those described for the ST12 clue above as 
items K. and L. 

If the Reinforcing Clue is ISG3 (unequal face sheet 
sandwich reinforced by a general angle isogrid 
construction)'; 

A. through E. The items contained on these cards 
are those described for the ISGl clue above. 

F. through H. The items contained on these cards 
are those described, for the ST13 clue above 
as items K. through M. 

If the Reinforcing Clue is STIC the following Hooke’s 
Laws will be used by the program for the description 
of the shell wall: 


k sets of 


5 sets of 
2 sets of 


5 sets of 
3 sets of 




- “lA - “to 

o 

^o 



+ ^ 

- C - N 

22 

12 0 

o 

22 0 T4> 


N = N = K y 

^33 M 
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5E14.7 

5E14.7 

5E14.7 

5E14.7 

5E14.7 

'5E14.7 

5E14.7 
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^0 ^ 11^0 ‘^ 11^6 ■ ^0 

o 

^<t> ^ ^22^4) '*’ ^12^0 ^22 ^4> " ^(j> 


M. 


"^ 04 . "^®33^4>e 


Therefore the input is (see Ref. 4 Appendix A); 


A. 


C. 


E. 


F. 


G. 


H. 


I. 


J, Initial, intermediate and fin 

at points defined by table of (|) or s values' 

Note ; No plasticity is allowed for this option. 

If the Reinforcing Clue is RWAF the following Hooke’s 
Laws will be used by the program for the description 
of the shell wall; 

^11 ^12 ^4> ■ ^11 ^0 " ^15 ^4) ' ^T0 

0 0 Y 


values of 
or s values. 

5E14.7 

values of 
or s values. 

5E14.7 

values^ of 
or s values! 

5E14.7 

values of 
or s values; 

5E14.7 

values of 
or s values. 

' sign convention. ) 

5E14.7 

values of 
or s values. 

' sign convention.) 

5E14.7 

values of 
or s values! 

* sign convention.) 

5E14.7 

values of 
or s values. 

5E14.7 

values of 
or s values. 

5E14.7 

values of 
or s values! 

5E14.7 


*^4> ^2 ^4)^ ^12 ^0^ ■ ‘^15 ^0 " ^22 \ " ^T4> 


%e ^33 ^4 


^^16 ^4)0 
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o 
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“ ^11 ® i 2 ^11 ^6 *^15 ^ 4 > ■ ^0 

^ O O 

“ ^22 ^(|i ^12 ^0 ^15 ^ 0 ^ °22 ^ 4 )^ ■ ^4 

%e “ ‘^^33 ^4>e ^i6 ^4»e 

O 

Therefore the input is (see Ref. 4 Appendix A): 

A. through J. The items contained on these cards are 10 sets of 
those described for the STIO clue above. 

K. Initial, intermediate and final values of C^(- e,t 
points defined by table of 4> or s values. 

, L. Initial, intermediate and final values of at 
points defined by table of 4 or s values. 

Rote ; No plasticity is allowed for this option. 

7. Table of Reinforcement Temperatures 

This card set is included only if the Reinforcement 
Clue (Item D page 2-39) is STll, ST12, ST13, RNAl, 

RWA2, RWA3, ISGl, ISG2, or ISG3 and the Thermal Clue 
(item E page 2-4o) is TEST or THGn. The input is 
presented below as a function of the Reinforcement 
Clue. 

If the Reinforcement Clue is STll, ST12, or ST13: 

A. Initial, intermediate and final values of the 
stringer temperature, at points defined by 
table of ^ or s values. 

B. Initial, intermediate and final values of the 
ring (smeared) temperatixre, at points defined 
by table of ^ or s values. 

If the Reinforcement Clue is RWAl, RWA2, EWA3, ISGl, 

ISG2, or ISG3: 

A. Initial, intermediate and final values of the 
reinforcement (waffle or isogrid) temperature, 
at points defined by table of <P or s values. 

8. Table of Mass Densities (Dimensions = weight /volume/ 
acceleration due to gravity. ) 


5E14.7 

5E14.7 

!?E14.7 


5E14.7 

5E14.7 

5E14.7 
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This card set is included only if the Shell Rotation 
Clue (Item K page 2-l8 ) is set to unity, i.e, if 

the shell undergoes rotation. In such a case, the 
contents of these cards (pp to 6 cards per item below) 
are dependent upon the clues registered on the Master 
Clue Card. If the shell to be described contains no 
reinforcing, the pertinent clue is item 4c, the sheet 
clue. The input is presented below as a function 
of the Sheet Clue. 

If the Sheet Clue is SING (single sheet construction): 

A. Initial, intermediate and final values of the wall 5®l4.7 

mass density at points defined by table of cp or s 
values. 

If the Sheet Clue is EQUA or UREQ (sandwich 
construction): 

A. Initial, . intermediate and final values of the 5^14.7 

core mass density at points defined by table of 

cp or s values. 

B. Initial, intermediate and final values of the face 5^14. 7 

sheet mass density at points defined by table of 

cp or s values. 

If the shell is reinforced, the Sheet Clue will be 
BIAN, In this case it is the following, or 
Reinforcement Clue (item 4 d) which will determine 
the contents of card series 8. ^he reinforced shell 


input is presented below as a function of the 
Reinforcement Clue. 

If the Reinforcement Clue is STll (single sheet 
reinforced by rings and/or stringers): 


A. Initial, intermediate and final values of the 
base shell mass density at points defined by 
table of cp or s values. 

5E14.T 

B. Initial, intermediate and final values of the 
stringer mass density at points defined by 
table of cp or s values. 

5E14.7 

C. Initial, intermediate and final values of the 
ring mass density at points defined by table 
of cp or s values. 

If the Reinforcement Clue is ST12 or ST13 ^sandwich 
reinforced by rings and/or stringers): 

5E14.7 

A. Initial, intermediate and final values of the 
base shell core mass density at points defined by 
table of cp or s values. 

5E14.7 
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B. Initial, intermediate and final values of the 
base shell face sheet mass density at points 
defined by table of cp or s values. 

C. and D. The items contained on these cards are 2 sets of 

those described for the STll clue as items 

B. and C, 

If the Reinforcement Clue is EWAl or ISGl (single 
sheet reinforced by a rotated waffle or isogrid): 

A. Initial, intermediate and final values of the base 
shell mass density at points defiied by table of 

cp or s values. 

B. Initial, intermediate and final values of the 
reinforcement mass density at points defined by 
table of 9 or s values. 

If the Reinforcement Clue is RWA2, RWA3, ISG2 or ISG3 
(sandwich reinforced by a rotated waffle or isogrid): 

A. Initial, intermediate and final values of the base 
shell core mass density at points defined by table 
of cp or s values. 

B. Initial, intermediate and final values of the base 
shell face sheet mass density at points defined by 
table of cp or s values. 

C. Initial, intermediate and final values of the 
reinforcement wass density at points defined by 
table of cp or s values. 

If the Reinforcing Clue is STIC or RWAF the program 
does not have any crossection description available 
to it for calculations. Thus no plasticity analysis 
is allowed, and the required input in card set 8 for 
these clues is (see Refs. 1, If and 5): 

A. The meridional distributed lead on the crossection 
due to rotation at points defined by table of cp or 
s values. For example for a layered section this 
can be expressed as: 


9 


^t^ cos cp 

i=l 


B. The normal distributed load on the crossection due 
to rotation at points defined by table of cp or s 
values. For example for a layered section this 
can be expressed as: 



i=l 


5E1^.7 

5Ellf.7 

5E14.7 

5E14.7 

5E14.7 

5Ellf.7 

5E14.7 

5E14.7 

5E14.7 
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C. The circumferential distributed moment load on 
the crossection due to rotation at points 
defined by table of cpor s values. For exan^le 
for a layered section this can be expressed as: 


me = 


§ 




'i"^i-i«'iC0S9 


Column 


Note: The sign convention for the above loads is 

as shown in Figures 2-11. 

9. Loading Clue Card 

The contents of this card are numerical clues which 
alert the program to the types of loads that exist 
on the segment. In addition, if the clue indicates 
that some load does not exist, the appropriate 
cards in series 10 which would ordinarily contain 
the numerical values of this load are omitted from 
the sequence. In the case of a rotating shell the 
inertia loads are calculated internally by the program. 


The appropriate clues are as follows: 

A . Thermal Clue 

If there are no thermal loads (Item Ue is 
NOTH) the clue number is zero (O). 

If there is a standard thermal variation through 
the thickness (item UE is THST) the clue number 
is four (U). 

If the temperature is constant through the 
thickness (Item UE is THCN) or if the in- 
homogeneous option is used (Item UE is THIN) 
the clue number is one (l). 

B. Circtunferential Load Clue (f ) 

0 

If there are no circumferential loads, then the 
clue number is zero (O). 

If there are circumferential loads, then the clue 
number is one (l). 


Format 

5E1U.T 


II 


II 
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C. Meridional Load Clue (f. ) 3 

If there are ^ meridional loads , then the 
clue number is zero (O). 

If there are meridional loads , then the clue 
number is one (l). 

D. Normal Load Clue (f^) U 

If there are rxo normal loads , then the clue 
number is zero (O), 

If there are normal loads , then the clue 
nximber is one (!)• 

E. Circumferential Moment Load Clue (m.) 5 

6 

If there are no circumferential moment loads, 
then the clue number is zero (O). 

If there are circumferential moment loads , then 
the clue number is one (l). 

F. Meridicmal tfoment due (m ) 6 

CP 

If there are no meridional moment loads , then 
the clue nuit±)er is zero (O), 

If there are meridional moment loads, then the 
clue number is one (l). 

G. Any alphameric information (load description) 7-70 

10. Table of Applied Loads (see Figures 2-lla, b for 
sign convention). 

The appropriate card sequence is given below as a 
function of the Loading Clues on card 9- If ^^e Thermal 
Clue is one (l): 

A. Initial, intermediate and final values of the 
temperature of the shell at points defined by 
table of (|i or s values. (These values will be 
used either for a thermal problem where there 
is no thermal variation through the thickness 
{Clue = THCN}, or to calculate varying material 
properties along the shell for an inhomogeneous 
problem { Clue = THIN } . ) 


II 


II 


II 


II 


16A4 


5E1U.7 
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If the Thermal Clue is four (U): 

A. Initial, intermediate and final values of the 
temperature at points defined by table of 
or s values. (The subscripts "nm" indicate 
temperature location - see below.) 


h. 

1 



B. Initial, intermediate and final values of the 
temperature T. at points defined by table of 

or s values; 

C. Initial, intermediate and final values of the 
temperature T at points defined by table of 

or s values • 

D. Initial, intermediate and final values of the 
temperature T at points defined by table of 
()) or s values . 

If the Thermal Clue is zero (O), the above cards are 
omitted. All temperature values above refer to the 
base shell. 

If the Circumferential Load Clue is one (l): 

E. Initial, intermediate and final values 
of the circumferential loads f^ at points 

defined by table of cp or s values. 

If the Circumferential Load Clue is zero (O), cards 
E are omitted. 


5E1U.T 


5E1U.7 

5E1U.T 

5E1U.7 

5E14.7 


/ 


2-^3 







I 


Figure 2-llc. Shell Element Geometry and Displacements 
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If the Meridional Load Clue is one (l): 


F, Initial, intermediate and final values of the 
meridional loads f at points defined by table 
of (|) or s values. ^ 


If the Meridional Load Clue is zero (O), cards Fare 
omitted. 


If the Normal Load Clue is one (l): 

G. Initial, intermediate and final values of the 
nonnal loads f at points defined by table of 
(|) or s values . ^ 

If the Normal Load Clue is zero (O), cards g are omitted. 


If the Circumferential Moment Load. Clue is one (l): 

H. Initial, intermediate and final values of the 
circumferential moment loads m^ at points 
defined by table of $ or s values. 


5E1U.7 


5E1U.T 


5E1U.T 


If the Circumferential Moment Load Clue is zero (O), 
cards H are omitted. 

If the Meridional Moment Load Clue is one (l): 

I. Initial, intermediate and final values of 5E14.7 

the meridional moment load m at points 

' CD 

defined by table of cp or s values. 

If the Meridional Moment Load Clue is zero (O), 
cards I are omitted. 

11. Segment Topology Cards 


A. 

Segment number 

Number of the segment under consideration. 

1-5 

15 

B. 

Joint (i) 

6-10 

15 


Joint associated with i'*^^ end of the segment 
(TIC). 



C. 

Joint (j) 

11-15 

15 


Joint associated with the end of ^the 
segment (STOP). 

Since within a region the segments are all 
singly connected, the segment J oint numbers 
should be in adjacent numerical pairs. That 
is, if joint (j) is 6, joint (i) could only 
be 5 or 7- This is true only within a region. 


2«57 



Column 


Format 


In addition, the initial joint of each 
region must he 1 in segment topology 
numbering, and, the final joint of each 
region must be the last (highest) number 
in the segnent topology numbering (see 
Figure 2-12). The coordinate or s 
increases from TIC to STOP, i to j. 

The user is again advised to see 
Figures 2-3 to 2-9. 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus 

(See General Notes - Data Debugging) 

INTRA-REGION DISCRETE RING CARDS 

These cards, if any exist (region introductory card, 
item IC), are placed at the end of all the segment 
data for the region. They contain the following 
information for each ring (in groups of 5 cards): 

1. Ring Property Card 

A. Segment joint number to which the ring is 1-2 

attached. If there is also a- radial 
discontinuity at this location necessitating 
a kinematic link, the joint number of the 
ring is the independent joint number of the 
link (see Figure 2-13a). 


B. Ring extensional stiffness (EA) 3_X6 

C. Ring bending stiffness about centroidal y 17-30 

axis (Ely). See Figure 2-lk. 

D. Ring cross-bending stiffness about centroidal 31-if4 

axes (EI ) . 
xy' 

E. Ring bending stiffness about centroidal 59"'^2 

X axis (EIjj.) • 

F. Ring Young's Modulus 45-58 

2. Ring Plasticity Card 

A. Ring Geometry Description Clue 1-4 


A specific clue word is necessary to describe 
the ring shape. The possible clue words and 
the associated shapes are presented in 
Figure 2-15. 


in l -80 


12 

E14.7 

E14.7 

E14.7 

Ei4.7 

ei4.7 

A4 
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z. 


z 



No Links Necessary 



Independent Joint 
of kinematic link 
and therefore Ring Joint 
number . 

Kinematic Link Necessary 


(a) Segment or Simple Region Rings 


Z 


kinematic link 1 



Independent Joint for all 
kinematic links and 
therefore Ring Joint 
number. 


(b) Possible Region Multi-connected Rings 


Figure 2-13. Discrete Ring Topology 
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I. 





Figtire 2-l4. Discrete Ring Geometry 
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Ring Orientation Angle 7-20 ElU^7 

Angle (in radians ) at which ring is oriented, 

Romally this entry is the s^e as the "phi 
input" for a curved shell; or is equal to 
the angle input for a cone (see Figure 2-8, but 
cp range here is between 0 and 3tt radians [no 
negative angles in inverted cone]); or is equal 
to tt/2 for a cylinder. 

Notes ; 1, In order to "invert” the ring, i,e., 
tvorn it upside down from the position shown in 
Figure 2-15^ as it rotates clockwise, and attach 
it to the inside part of the shell, add n to the 
angle above , 

2, If the ring is oriented at another 



angle rather than tangent to the shell wall, 
the angle should be entered accordingly above. 



c. 

Ring Yield Stress 


21-3k 

E14,7 

D. 

Ramberg •» Osgood constant 
ring (see Ref. 1). 

"s" for the discrete 

35 >48 

E14.7 

E. 

Ramberg - Osgood constant 
ring (see Ref. l). 

"n" for the discrete 

49 *"62 

E14.7 


3. Ring Location Card (see Figure 2-l4) 


A. 

Ring thermal coefficient of expansion (o^)» 

1-12 

EI 2.5 

B. 

Ring centroidal radius (r^). 

13-24 

EI 2.5 

C. 

Distance (±) between ring centroid and shear 

center in the x direction (x ). 

c 

25-36 

E 12.5 

D. 

Distance (±) between ring centroid and shear 
center in the y direction (y^)- 

37-48 

E 12.5 

E. 

Offset (±) of ring shear center from attached _ 
shell Joint (see item lA) in the x direction (x). 

49-60 

E 12.5 

F. Offset (±) of ring shear center from attached 

shell Joint (see item lA) in the y direction (y). 
Note; The shear center fS ) is the origin of the 
X, y axes shewn in Figure 2-l4, 

Ring Thermal Description Card (See Figure 2~l4) 

61-72 

E 12.5 
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A. Innermost point temperature (T^) 

l-l4 

E14.7 

B. Outermost point ten^jerature (T^) 

15-28 

E14.7 

C. Radius of revolution at the segment joint which 
is used to identify the ring (item lA). 

29-42 

E14.7 

D. Ring mass density (P ) 

43-56 

E14.7 

E. Kferdening Law Clue 

61-64 

A4 


Different hardening laws can be used 
for each ring (as well as each segment). If 
isotropic hardening is to be used the clue 
is ISOT . If kinematic hardening is to be used 
the clue is KUTE . If perfect plasticity is 
to be used the clue is PERF . For this last 
clue word items 2D, E( above) are set to zero . 

5. Ring Geometry Card (see Figure 2-15) 

The ring geometry is specified on this card (as many 
items as necessary) as a function of the Ring Geometi^ 
Description Clue (itm 2A). 


If the clue word is ZSEC, ISEC, or CHSC: 


A. 

Length dimension, a 



1-12 

E 12.5 

B. 

Length dimension, b 



13-24 

EI 2.5 

C. 

Length dimension, c 



25-36 

EI 2.5 

D. 

Thickness dimension, t^ 



37-48 

EI 2.5 

E. 

Thickness dimension, t^ 



49-60 

E 12.5 

F. 

Thickness dimension, t^ 



61-72 

EI 2.5 

If the clue word is HTRI, CRRC, TSEC, or 

RASC; 



A* 5 

B. These items are the 
above for the ZSEC 

same as those identified 
clue . 

1-24 

2E12.5 

C., 

D. These items are those identified 

above for 

25-48 

2E12.5 


the ZSEC clue as items D and E. 
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001x31011 Format 

If the clue word is HEIEC; 


A-D. These items are the same as those identified above 
for the HTRI clue 

1-48 

4e12.5 

E. Thickness dimension, t^ 

49-60 

EI2.5 

If the clue word is SKEC: 



,A. Length dimension, b 

1-12 

EI2.5 

B. Length dimension, a 

13-24 

E12.5 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 

mimm in 

H 

CD 

0 


(See General Notes - Data Debugging) 
INTRA-REGION KINEMATIC LINK CARDS 

These cards> if any exist (region introductory card, 
item 1B), are placed at the end of all the discrete 
ring data for the region. They contain the following 
inf omation : 


A. Joint (j) dependent joint 

1-2 

12 

B. Joint (i) independent joint 

For intra-region kinematic links these joints 
must be in consecutive decending order. Thait 
is, joint (j) should always be greater than 
joint (i) by one. 

3-4 

12 

C. Angle 7 in radians (see Figure 2-16) 

7 cannot equal 0 or tt. 

The angle 7 describes the orientation of the 
link; it is the inclination angle of the link 
from the vertical (Z axis). 

The nxomber of kinematic link cards must 
equal the number specified in item IB 9 
of the region introductory card. 

5-19 

E14,7 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 
(See General Notes - Data Debugging) 

minus 

in 1-80 


REGION JOINT CONTROL DATA 

These cards are placed at the end of all the data for 
all regions. 

1. Joint Control Data Card 
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Figure 2-16. Kinematic Link 



Colmn 


Format 


A. Number of region joints 

Total number of region joints (Max* - 20)* 

B* Number of discrete rings 

Total number of discrete rings between 

regions. 

C. Number of kinematic links 

Total number of kinematic links between 

regions. 

2. Discrete Ring Cards (inter-region, if any exist) 

A. Ring Property Card 

a. Region joint number to which the ring is 
attached. If there are also radial 
discontinuities at this location neces- 
sitating kinematic links, or multi- 
connections, the joint number of the ring 
is the independent joint number of all 
the links (see Figure 2 -13b). 

b. Ring extensional stiffness (EA) 

c. Ring bending stiffness about centroidal 
y axis (EI^). See Figure 2-l4. 

d. Ring cross-bending stiffness about 

centroidal axes (El ). 

^ xy^ 

e. Ring bending stiffness about centroidal x 
axis (EI^). 

f. Ring Young ^s Modulus 

B. Ring Plasticity Card 

a. Ring Geometry Description Clue 

A specific clue word is necessary to describe 
the ring shape. The possible clue words and 
the associated shapes are presented in Figure 
2 « 15 . 


l>5 

15 

6-10 

15 

11-15 

15 


1-2 

12 

3-16 

E14.7 

17-30 

E14.7 

31-UU 

E14.7 

45-58 

E14.7 

59-72 

E14.7 

1-4 

a4 
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Column 


Format 


b« Ring Orientation Angle 7-20 El4,7 

Angle (in radians ) at which ring is oriented* 

Normally this entry is the same as the "phi 
input" for a curved shell; or is equal to 
the angle input for a cone (see Figure 2-83 
but cp range here is between 0 and Stt radians 
[no negative angles in inverted cone]); or 
is equal to tt/2 for a cylinder.* 

Notes: 1* In order to "invert" the ring, 

i*e*, turn it upside down from the position 
shown in Figure 2«15j as it rotates clockwise, 
and attach it to the inside part of the shell, 
add tt to the angle above* 


2. If the ring is oriented at another 
angle rather than tangent to the shell wall, 
the angle should be entered accordingly above. 


c. Ring Yield Stress 

21-34 

E14.7 

d. Ramberg - Osgood constant "s" for the discrete 
ring (see Ref. l). 

35^48 

E14.7 

e. Ramberg - Osgood constant '^n" for the discrete 
ring (see Ref. l). 

49-62 

E14.7 

Ring Location Card (see Figure 2-l4) 



a. Ring thermal coefficient of expansion 

'“r)- 

1-12 

E12.5 

b» Ring centroidal radius (r ), 

c 

13-24 

EI2.5 

c* Distance (±) between ring centroid and 
shear center in the x direction (x ). 

25-36 

EI2.5 

d« Distance (±) between ring centroid and 
shear center in the y direction (y^)« 

37-48 

EI2.5 

e. Offset (±) of ring shear center from 

attached shell Joint (see item Aa) in the 
X direction (x). 

49-60 

E 12.5 

f. Offset (±) of ring shear center from 

61-72 

EI2.5 


shell Joint (see 'item Aa) in the y 
direction (y*). 

Note : The shear center (S) is the origin 

of the x,y axes shown in Figure 2-l4. 



Column 


Format 


D, Ring Thermal Description Card (See Figure 2-l4) 


a. 

Innermost point temperature (T^). 

l-l4 

ei4.7 

b. 

Outermost point temperature (Tq). 

15-28 

ei4.7 

c* 

Radius of revolution at the region Joint 
which is used to identify the ring (item 

24-42 

E14.7 


Aa). 



d. 

Ring mass density (p ) . 

43-56 

E14.7 

e. 

Hardening Law Clue 

61-64 

a4 


Different hardening laws can be used for 

each ring (as well as each segment). If 

isotropic hardening is to be used the clue 

is ISOT . If kinematic hardening is to be 

used the clue is KINE . If perfect plasticity 

is to be used the clue is PE^. For this 

last clue word items Bd, e (above) are set to zero. 

E. Ring Geometry Card (see Figure 2-15) 

The ring geometry is specified on this C8u*d 
(as many items as necessary) as a function 
of the Ring Geometry Description Clue (item 
Aa). 

If the clue word is ZSEC, ISEC, or CHSC: 


a. Length dimension, a 1-12 E12.5 

b. Length dimension, b 13-24 E12.5 

c. Length dimension, c 25-36 E12.5 

d. Thickness dimension, t^ 37-48 E12.5 

e. Thickness dimension, tg 49-60 E12.5 

f. Thickness dimension, t^ 6l-72 E12.5 

If the clue word is HTRI, CRRC, TSEC^ or RASC: 

a. , b. These items are the same as those 1-24 2E12.5 

Identified above for the ZSEC clue. 

c., d. These items are those identified above 25*^48 2E12.5 

for the ZSEC clue as items d and e. 
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Cnlunn 


Format 


If cXu£ word 1 e HEiSC : 


a-d* These items are ttw sanse as those identified 
above for the EEBI clue 

iJtB 

4ELE.5 

0. 

Tttickaess dimension, 



If 

tba clue is EEtSC: 



It. 

Length dlmensltmj b 

1-12 

£12.5 

b. 

Length dimension 1 a 

13^^ 

E12.5 


D^-5-H S-£-P^-B-IV-T-0-R C^-B-D mlnuB In 1-00 

[3*e General Wo tea - Data Debugging} 


/ 


/ 


3, (tlnfrnatlc L1 j])£ Cards ( intar-reftion* if any exist) 


Joint {j) dflponiBnt joint 

1-2 

12 

B- Joint fi) indepenjiant ;Jo±nt 

3-^ 

12 


Far kinenatic links between regiaaa ttvero 
are no rOStrletlOnO upon Joint Kunnbering. 


The only restriction is ttlat betwoon succOEElve 
kinematlo link data oardE the (j) ^cltlt entry 
should be In incrcaaiug order (not □eoessarlly 
ooneocutivo). For oxanqplej 


Host he 
in 

IncreaEllEig 

order. 


i i 

3 a 
h 2 
9 1 
12 6 


y 


C. Angle r In radians (see Figure 2-^16) 


T CMiigt e^ual 0 or 

_The— g^Le y describes the orientation of ttM 
linkj it is the inclination anSl^ of the 
link from the vertical (S axis). The tiUmber 
of klKianatle link oerde must equal the number 
specified in item 1C of the Joint Control 
Data Card. 


5-19 
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Colmm 


Fprn^t 


D-A-5-H S-B-P-A-R-A-T-O-K C-A-H-JD nlntie In 1-30 

(Bee General Hates - D^ta Dst^jLgiatnfi) 

Boundary Condition Carda {Joint Slftta - one card 
]?<*r region .1oint) 


Joinrt Huiiiher 

1-2 

la 

Joint component e&rMsJitiona on: 



i) Ap 


P£.0 

e) ^ or ffiee Ri^urc £-17? 

5-6 

P2.0 

3 ) 

7-8 

P2.0 

4) !1, 

9-10 

F2^0 

Anfil^i a Ln radians 

11-2^ 

El^^.l 


Co be Used only in eonj Uiiction with & S or 
3 code. 

Note : rtiETe nnat lie as cany boimdary con- 

dition oardB an there are Joints as indicated 
in Item lA of tho Joint Control Etata Cardi 


There are ^ different codes that are used 
to proiserlbe Joint component conditions + 
They lure: 

a- 0 ■ no diaplitcem^nt allowed. 

b- 1 - displacement alloved In the 
indicated direction. 

C. 2 ^ ^ an'l ^ we rotated throd^ ad 
angle Qf"(T/& - ■e) *nd becewe i- 
ani ^ reapoctlvcly, vtilLc a 
ilEplftcement ie allOiVfid lu the 
direct 
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Matrix 


Rotation Code 


Matrix 


Rotation Code 


# 

Matrix 


f 


+1 1 

0 

1 

1 

0 



0 j 

sin a 

1 

1 

cos a 


> * 

0 I 

-cos a 

1 

1 

sin a 

l“ej 


0 1 

0 

1 

1 

0 


1 

2 0 


1 

✓ N 


+1 1 

0 

1 

1 

0 

"z 

c 


0 { 

sin a 

1 

0 

> * 

0 

-cos a 

1 

1 

0 



0 

0 

1 

1 

+1 









^7. 


.“e. 

> “ 


+1 

0 

0 

0 


0 3 

0 

cos a 
sin a 
0 


0 

0 

0 

+1 



Figure 2-1713. Provision for Local Rotations 
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d. 


3 = 


\ and Aj^are rotated through an angle 
or (ir/2 - a) and become A|^ and ZV 
respectively, while a displacement 
is allowed in the direction. 


See Figure 2-17 for a geometric 
explanation of codes 2 and 3. 


When using rotation codes: 

Code 2 can exist only as Z^ coding. 

Code 3 can exist only as Z^ coding. 

Codes 0 and 1 can appear in either 
column 4 or column 6, in addition 
to columns 8 and 10. Thus, there 
are twelve possible boundary con- 
ditions when rotation codes are 
used, (a - (p for table below.) 



Free edge (possible 
to apply shear and/ 
or membrane loads) 

Z^ = 0, normal 
support (possible 
to apply membrane 
load) 

= 0, membrane 
support (possible to 
apply shear load) 

Z4p, % free 

1,2, 3,1 

1^ 2^0^ 1 

1,0, 3,1 

Qq fixed 

o,2,3,o 

0,2, 0,0 

0,0, 3,0 

Z^ fixed 
free 

0,2,3,1 

0,2, 0,1 

0,0,3,1 

free 

Q fixed 
0 

1,2,3,0 

1,2,0,0 

1,0, 3,0 
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Column 


Format 


Apex boundary conditions ; 

Since the closed apex ((ji * O) angle 
is not acceptable input, the apex 
boundary conditions must be simulated 
at a small ^ angle. These boundary 
conditions vsiry per Fourier harmonic, 
and are as follows: 


n 

Ap 



"e 

angle 

0.0 

0 

0 

3 

0 

a - ^ 

1.0 

1 

2 

0 

1 

II 

>2.0 

0 

0 

6 

0 



General Notes ; 

1} To establish a datum for measuring 
displacement, free body motion must be 
eliminated from the structiare. This should 
be accomplished by suitably applied boundary 
conditions. 

2) The ability of a dependent Joint in a 
kinematic link to prescribe motion independ- 
ently should be removed by setting all 
boundary conditions of that Joint to zero. 
See pages 2-65 and 2- 8H 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 

(See General Notes - Data Debugging) minus in l-80 

JOINT IPAD DATA 

1. Load Control Data Card 

A. Number of Joint Loads 1-4. j4 

Total number of Joint loads in analysis. 

(Line loads can only be applied to region 
Joints. ) 

B. Any alphameric information (load description) 5-68 16 a4 

Rote ; If there are no Joint Loads for the 
structure, card 1 of the JOINT LOAD DATA 
is blank and card set 2 is omitted. 
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Column 


Format 


2. Joint Load Cards (as many as in item LA above) 


Enter unity in colTmin 5 

1-5 

15 

Row Identification 

6-10 

15 


The identification is the location of the 
degree of freedom at which the load is 
applied. This is obtained _^by counting the 
non-zero codes entered in the Boundary 
Conditions Cards, starting with Joint 1; 

T, Z, R, fig. Joint 2; T, Z, etc., 

and stopping at the joint and degree of 
freedom where the line load is to be 
applied. The location number of this 
degree of freedom is the information 
necessary. 

C. Applied Joint Load 11-24 E14.7 


The input is 2irr times the running load 
in lb. /in. In tSe particular case of the 
axial axisymmetric load, this is simply 
the net force. For sign convention see 
Figures 2-11, 2-1 7 and 2-1 d 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in l-80 

(See General Notes - Data Debugging) 


LOAD CYCLING CARDS 

If only a single loading is being investigated these cards are 
omitted. For changes of load distribution, magnitude, or mat- 
erial property tables with time, the appropriate cards below 
are included. (See also item 2H in the General Introductory 
Cards p 2-i8 ). These cards are repeated to amoimt described 

in Item 2E on p 2-l6 

1. Load Cycling Program Control Card 

A. Intermediate print clue 8-9 12 

(See item 2D General Introductory Cards p 2- l6). 
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Column 


Format 


B. dumber of load steps 12«17 

Number of loading increments to be made for the 
next loading cycle. If the loading reverses 
in direction^ this number is input as being of 
oppos ite sign to the same entry for the previous 
loading cycle. 

C. Load step print increment for next loading cycle* l8«*2l 

( See Item 2G General Introductory Cards p 2-l8 , ) 

D* Load amplitude change clue 24-25 

If the previous loading cycle load amplitudes are 
to be reused^ and the only changes involved are 
load directions or step sizes (controlled by item 
C above) or material property tables ^ the clue is 
zero. If the load amplitudes are also changed^ the 
clue is five , and new load cards will be expected 
to be read. It should be noted that load amplitudes 
are sequential. For further clarification see 
Figures below. 


F6.0 


f4.o 

12 
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Column 


Format 


E* Material Table Change Clue 26-27 12 

If the previous loading cycle material property 
tables are to be reused the clue is zero . If 
new material property tables are supplied the 
clue is one. 

Note ; At the first load reversal in direction 
the material property tables must be re-read 
with all the Ramberg - Osgood ”s” constants 
doubled in value. This is done only 
at the first load r ever sal and the new tables 
can be reused thereafter. If the constants 
above are all zero this can be disregarded. 

This is true for all reinforcement as well as 
the shell. 

F. Angluar velocity in rad. /sec. for next loading 59“72 El4.7 

cycle. 

Since the shell rotation clue (item 2K in the 
General Introductory Cards p 2-18 ) is constant 

for the total run^, the above item is significant 
only if that clue was initially set to one. In 
order to apply or delete an angular velocity in 
the middle of a loading history^ the velocity 
values for other load cycles may be set to zero. 

2. Material Property Tables for Next Loading Cycle 

If the preceeding Material Table change clue was set 
to one a complete new set of material property tables 
must be provided. Note the special item above for the 
Ramberg - Osgood "s” parameter. The names of the 
material tables may not be changed. The necessary 
formats and card sequences are provided on p 2-21 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-B minus in l-80 

(see General Notes - Bata Bubugging) ■ 
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Column 


Format 


If the preceedins Material Table Change Clue was set to zero. 
item. 2 above and the separator card are omitted. 

3. Shell loading data for next loading cycle 

If the preceeding Load Amplitude Change Clue was set 
to five , the following card sequence must follow. 

A . Segment Loading Clue Card for next loading cycle 

a. Load clues for temperature, circumferential, 1-6 6ll 

meridional, normal, circumferential moment, 
and meridional moment load, for next loading, 
cycle . 

Note ; Variations in types of loadings as 
well as magnitudes are allowed per cycle. 

The thermal clue, however, is fixed and 
must remain constant. Thus AT = 0 may have 
to be input as appropriate. 

B. Table of Applied Loads for next load cycle (see 5E14.7 

segment card set 10). These cards contain the 

next set of loads as described previously in the 
segment data according to the Load Clues above. 

Hote ; The above data sets (3A and 3B) are re- 
peated until all the segments in the first 
region are con^jleted for the next load cycle. 

C. Region Ring Load Cards for next loading cycle. 


a* Ring Yield Stress 

l-l4 

E14.7 

b. Ring Ramberg - Osgood value (see note 

about doubling this value for first load 
reversal, above.) 

15-28 

E14.7 

c. Ring Ramberg - Osgood ”n" value. 

29-42 

E14.7 

d. Innermost point temperature (T^). 

43-56 

E14.7 

e. Outermost point temperature (T^). 

56-70 

E14.7 


ITote : The above data set (3C) is repeated 

until all the rings in region one are completed 
for the next load cycle. 


The above data sets (3A through 3C) are then 
repeated until all the segments and rings in 
all the regions are completed for the next 
load cycle. 
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/ 

Coltann Fonnat 

The above data set (3C) is then repeated until 
all the structure rings are completed for the 
next load cycle. 

D. Joint Load Control Data Card for next loading cycle. 

a. Number of Joint Loads 1-4 l4 

b. Any alphameric information 5-68 16 a4 

Note; If no joint load changes exist in this load 
cycle, the above card is blank and card set 3E (below) 
is omitted. 

E. Joint Load Cards for next loading cycle 


a. 

Enter unity is column 5 

1-5 

15 

b. 

Row identification (see previous Joint Load 
Card description) 

6-10 

15 

Note: Locations of line loads may change from 

load cycle to cycle. 



c. 

Applied Joint Load 

11-24 

E14.7 


The input is 2TTr^ times the running load in 
lb. /in. for this load cycle (see Figure 2-l8). 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in l-80 

(see General Notes - Data Debugging) 

If the preceedlng Load Amplitude change clue was set to zero. 
item 3 above and the separator card are omitted. 

The above cycling cards (items 1-3) are repeated until all the 
load cycles input in item 2E on p 2-l6 are completed. 
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z 



Figure 2-1 9- Y-Joint (Distorted Geometry) 
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It is hoped that the user is now able to use the SIARS-2P program to good 
advantage. It is a powerful tool which will increase in value to the user as 
he uses it. One of the more complex areas of usage is the description of 
topology, especially when involved with rotation codes and Joint loads. An 
illustrative example of a Y Joint representation is therefore presented below. 
(See Figures 2-19 and 2-20 for the structure and idealization. ) Ihe idealized 
stmctufe contains four regions and two kinsaatic links. The Joints are numbered 
from 1 to 7. Membrane loading is applied to Joints 1 and 5 and the structure 
is supported by membrane action at Joint 7. All regions must be coupled. 

The second card in each region description (topology card) is as follows; 



The Joint Control Data card would contain a 7 in column 5 and a 2 in column 15. 

In this example, the restraints at Joints 1, 5 and 7 must be rotated from the 
fixed (global system) to a local system such that membrane action may be 
applied. In addition. Joints 3 and 4, and 6 and 3 are to be coupled with 
kinematic links. Thus, the motions of Joints 4 and 6 are dependent upon the 
motion of Joint 3. This dependence will be insured by using 2 kinematic link 
cards and setting the independent degrees of freedom of Joints 4 and 6 equal to 
zero. It should be noted in this particular case that the motion of Joints 4 
and 6 is not being equated to zero, but rather, the ability to prescribe 
motion independently is being removed. The required data has the following 
appearance. 

A. Kinanatic Link Cards (2 cards) 

Note * In a double link of the type 
shown in Figure 2-20, one Joint 
must be consistently independent 
(Joint 3 in example). 
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B. Bountary Condition Cards for a coupled torsional case (7 cards) 
Joint A ^ Angle 

1 1 2 3 1 

2 1111 

3 11 1 1 

4 0 0 0 0 

5 12 3 1 02 

6 0 0 0 0 

7 0 0 3 1 *^3 

The Load Control Data card would contain a 2 in column 4. 

The external membrane loads are applied to the structxire through the 
Joint Load Cards which, in this example, would appear as (2 cards): 

Problem No. Row Load 

1 2 ^ Membrane Load (lb. /in.) 

1 1^ ^ Membrane Load (lb. /in.) 
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SECTION 3 


OUTPUT INFORMATION 


The output of the STARS-2P program is straightforward, however a 
description is in order since the user should learn the significance 
of the various checks that are provided. It is important to point 
out that the output of the program will include a print-out of the 
input data. This gives the user the opportunity to check whether 
or not the input data was correct. In the detailed description of 
the complete output which follows, the user should refer to the output 
of the problem in Section h as an example « 

The title page of the output contains all the data from the General 
Introductory Cards, prominently placed, and needs no comment. The 
next page of the output contains the first region Identification 
Card in the center. The following output is then presented for each 
segment in this region (in order of appearance): 

1. Contents of segment Identification Card. 

2. Contents of MAGIC Control Card 

3. Contents of Geometric Description Card (Cards) 

4. Contents of Master Clue Cart 

5. The material property table used for the segment 

6. Crossection description table 

7. Temperature load table (if any) 

8^ Distributed load tables 

9» Segment influence coefficients (MAGIC output) 

10. Segment stiffness matrix 

11. Stiffness matrix symmetry check 

12. Segment load matrices 

13. Radius of revolution at i^^ and ends of segment 
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Item llj the stiffness matrix symmetry cheeky is a check upon the 
validity of segment sizing and the accumulation of round-off error. 

For perfect symmetry to exist, it is necessary to have zeros above 
the main diagonal, and zeros or ones below the diagonal. The amount 
of error induced by improper sizing or round-off is related to the 
amount that the off-diagonal terms in the, lover triangle differ 
from unity. An attempt should be made to keep the upper limit on 
this difference at one percent (maximum number in lover triangle 
should be 0.1010... E Ol). 

As mentioned previously, items 1-13 are repeated for all the segments 
within region one. The radius of revolution at every joint should be 
checked at corresponding joints of adjacent segments to make sure 
that proper coupling has been specified. At this stage in the output, 
the topology of the segments within the region, and the description 
of the intra-region rings and kinematic links is presented. 

Next the region matrices are presented. Given in order are: the 

region stiffness matrix, the stiffness matrix symmetry check, and 
the region load matrices. Again the numerical round-off, evident 
in the symmetry check, should be kept to a maximum of one percent 
( 0.1010... E 01 in lover triangle). The output to this point, that 
is, sets of items 1-13^ segment topology and links, and the region 
matrices, are now repeated as a group for each region within the 
structure, "When this is completed, the region topology is presented. 

The next items to be provided by the output are the descriptions of 
the inter-region rings and kinematic links and the boundary conditions. 
At external points of the structure these are physical boundary 
conditions. At internal points they merely state the fact that no 
restraint exists and the joint in question is free to move. The last 
column in this set, gives the angle o', which is zero unless a rotation 
code is indicated . It is important to refer to Figure 2«17once more and 
point out that O' represents a rotation of the coordinate system. 
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There are a variety of errors that can be made in submitting input data^ 

The STARS-2P program is set so that as an error is detected ^ (input or 
in the matrix calculations) the program is stopped. Therefore 3 to avoid 
delays in getting an answer 3 the SATELLITE-IP program was created to 
debug input data^ and to check for other possbile inconsistencies. The use 
of this program is described in Reference 3 9 and all the STARS error messages 
are contained in the listings therein. 

After having corrected any errors 3 so that it is now possible for the pro- 
gram to run to completion^ the problem solution can be discussed. From 
the above point in the output 3 further output will depend upon the inter- 
mediate print clue. If this clue is set for full output 3 all stiffness 
matrices will be printed for all segments 3 regions 3 and struct-ures in all 
print cycles specified. If the print clue is set for abbreviated print 3 
only that noted below will be printed: 

1. The applied line loads. 

2. The region end deflections starting from region joint 1. 

3- Tabulation of final print (per print interval) as shown in Table 3-1 • 
Below this tabulation will appear the numerical results. In addition 3 
if the crossection has gone plastic the following results will be 
available ; 

a. RPLA; if this indicator is zero the layer is elastic. 

if this indicator is + the layer is plastic, 
if this indicator is - the layer is elastic 
but has previously been plastic . 

b. Total stresses through the crossection. 

c. Plastic strains through the crossection. 

In item 3 above the following special cases may arise: 

1, The I 32 element will be zero for cylinders or cones, 

2. No stresses will be calculated for the STIO or RWAF options, 

4, A restarted run will provide the initial print and start item 3 at 
the next print interval. 

This completes the program output. The user is reminded to check continuity 
of stress resultants across segment boundaries where applicable. 
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TABLE 3“1. Tabulated Final Output Information 


1,1 

2,1 

3.1 

4.1 

5.1 

6.1 

PHI (RAD. OR IN,) 
EPSILON THETA 
U 
V 
W 

OMEGA THETA 

CO or s at which calculations are made 

circumferential strain 
9 

circumferential displacement 
meridional displacement 
normal displacement 

circumferential rotational displacement 
(rad) 

1,2 

DEGREES 

cp, expressed in degrees 

2,2 

EPSILON PHI 

e 9 meridional strain 
CD 

3,2 

K PHI THETA. 

^cp9’ specific twist 

4,2 

J PHI STAB 

non-linear effective transverse shear 
cp 

stress resultant 

5,2 

T PHI THETA 

T effective in-plane shear stress 

cdB 

resultant 

6,2 

OMEGA PHI 

Cl) 9 meridional rotational displacement 

CD 

(rad) 

1,3 

R ZERO 

r , radius of revolution 
0 

2,3 

GAMMA FHI THETA 

Y ^9 shear strain 
coO 

3,3 

N THETA 

N ^5 circumferential force resultant 

4,3 

M THETA 

M .^9 circumferential bending moment resultant 

5,3 

SIGMA THETA IK 

Op^in? circmferential stress on inside fiber 

6,3 

SIGMA THETA OUT 

CT^out^ circumferential stress on outside 
fiber 

1,4 

BASE THICKNESS 

basic shell wall thickness 

2,4 

K PHI 

k 5 meridional curvature 

3,4 

N PHI 

T 

N 5 meridional force resultant 

CD 

M , meridional bending moment resultant 

CT . . meridional stress on Inside fiber 
^xn 

4,4 

M PHI 

5,4 

SIGMA PHI IN 

6,4 

SIGMA PHI OUT 

a .5 meridional stress on outside fiber 

cpouv 

J 
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TABLE 3-1. Tabulated Final Output Information (Continued) 

1,5 

STEP 

numerical integration step size 

2,5 

K THETA. 

k , circianferential cunrature 
9 

3,5 

N PHI THETA 

N in-plane shear stress resultant 

^,5 

M PHI THETA 

M' , twisting moment resultant 
(00 

5,5 

TAU PHI THETA IN 

T ' . , in-plane shear stress on inside 
fiber 

6,5 

TAU PEC THETA OUT 

T - in-plane shear stress on outside 

fiber 

1,6 

EPSILON THETA IN 

Extreme edge strain for base shell (inner) 

2,6 

EPSILON THETA OUT 

Extreme edge strain for base shell (outer) 

3,6 

EPSILON PEC IN 

Extrene edge strain for base shell (inner) 

4,6 

EPSILON PEC OUT 

Extreme edge strain, for base shell (outer) 

5,6 

GAMMA PEC THETA IN 

Extreme edge strain for base shell ( inner) ^ 

6,6 

GAMMA PHI THETA OUT 

Extreme edge strain, for base shell (outer) 
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SECTION 1)- 


EXAMPLES OF INPUT AND SOLUTIONS 


This section contains the full input for two example test problems, and 
the full output of the second test problem for the STAES-2P program. The first 
input problem was chosen simply to illustrate a load cycling run data set-up. 

The test problem is presented in Figure 4-1. The loading involves both 
distributed and line loads. The Intermediate print clue has been set for 
abbreviated output and graphic displays of some output variables have been 
requested. The necessary input is presented on the following pages. 






PLATE 

UNDER CYCLIC LOAD 

4 4 2 

3 120 20 2 

NLPL 

1 

ALUM 

ISO? 


10.05 

E 06 

.33 



8*411? g 03 

3*967 

4*0 £ 03 

ALUi ISO? 

i0«03 £ 06 

*33 

8*41? E 03 

3.967 

4.0 £ 08 


100 PLATE 
112 
21 20 

.002 1.0 £<*04 .0006 

0.0 

ISO? ALUl SING THIC 

.0020 .02 

.2613 .2615 

000100 

12.000 £ 03 12.000 £ 03 

112 

100 PLATE 
2 2 3 

21 20 

.01780 1.0 £«04 ,005933 

0.0 

ISO? ALUM SING THIC 

.0200 .1875 

.2615 .2615 

000100 

12.000 E 03 12.000 E 03 
112 


100 PLATE 
3 3 4 

21 20 



KINE 


2 



21130 


E-04 


0422S0 


0 


l.O 

6.0 

ISOT ALUM SING THIC NOTH KINE 

•1875 «6l 

•2615 *2615 

112 


10 0 PLATE 
6 4 5 

21 20 


1.0 
0.0 
ISOT 
• 61 
• 2615 

1.0 E-04 .2 

ALUM SING THIC 

2.61 
• 2615 

NOTH 

KINE 

0. 

2 

1 

1 2 




5 

0 0 





10 10 0 
2 0 111 

3 0 1 11 

4 0 111 

5 0 0 1 1 


1 

1 1 -.150797 


4 4 2 3 -190 20 5 I 


ALUM 

ISOT 

10.05 

m 

o 

.33 

1.6834 

o 

UJ 

3.967 

4«0 

E 03 

ALUl 

ISOT 

10.05 

m 

o 


• 33 



1.6834 E 04 
3.967 

4.0 E 08 


4-4 



OOOlOO 

I9«000 E 03 I9«000 E 03 
OOOlOO 

19.000 £ 03 19.000 E 03 


I 

1 1 «*.238762 


442 3 190 20 0 
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The second test problem Involves monotonic loading of a structure 
similar to that shown in Figure 4“1. The thickness of the plate for 
this problem Is 0ol286" and the loaded area radius Is 0,18675*'. The 
applied load Increases monotonical ly to 5,000 ps I • The input and out- 
put for the problem is presented on the following pages. The apex 
segment was treated as an elastic plug In this problem because of a 
local stress singularity at the apex which Is dependent upon computer 
accuracy. Further discussions of this and other problem results are 
aval l^ble in Ref. 1 . 


k-6 



CIRCULAR PLATE TPI 


4 4 2 

1 

NLPt 

1 

ALUM 

ISOT 

10.5 

£ 06 

• 33 

6.634 

£ 03 

3.661 

4.0 

E 03 

ALUl 

ISOT 

10.5 

E 06 

• 33 


6«634 E 03 
3«66i 

4«0 E 08 


100 PLATE 
112 
21 20 

<»002 l.O 

0.0 

I SOT ALUl SING 

•0020 «02 

•1263 •1283 

000100 

0«5000 E 04 0«5000 E 04 
112 


E“04 •0006 

THIC NOTH 


KINE 


100 PLATE 
2 2 3 

21 20 

.01780 l.O 

0.0 

ISOT ALUM SING 

•0200 .1875 

•1283 .1283 

000100 

0.5000 E 04 0.5000 E 04 
112 


E-04 .005933 

THIC NOTH 


KINE 


100 PLATE 
3 3 4 
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21 20 

*21130 l«0 E-04 *042250 

OaO 


ISOT 

*18T5 

*1283 

1 

ALUH SING 

«61 
• 1283 

1 2 

THIC 

NOTH 

10 0 

4 

21 20 

PLATE 

4 5 



l«0 
0 « 0 

i«0 E-04 

• 2 


ISOT 
« 61 

• 1283 

1 

ALUH SING 

2*61 
• 1283 

1 2 

THIC 

NOTH 

5 

0 0 




1 0 I 0 0 

2 0 111 

3 0 111 

4 0 111 

5 0 0 1 1 

1 

1 I 125664 


KINE 


KINE 





UnSYmhETrK* ORTHoTRoPtC* An^U^5*S «<TH CoMrEIhG OP AT hOST SHE^j, REMOnS 

■ - - “SrAKS^tP ■ 

AS OF ^MUY l« 1973 


NUHBER OF SE^HENTS « NOHBER Of REGtOtJS w *4 NUHBER Of HaTERJAE pROPgRTt TABLES USgO 2 WO* Of tVCLEE 

lAROEOEFEECTTON PLASTIC RROBLEH 


rt 


CIRCULAR PLATE TPl 


For jnfoRHation call v* svalbonas 

P* OBILVIE 

H«’ teviwg. 


Bo* 


— r 


575*7701 



REliah hun&iii I 


THMe AftE l> ItqpiEtiri AltD i tENMAT^C ,.||Hl IllIfU fH|* 





PlippfIXT VitD 


<iC4P4C 

*VWO 


frHiaafl' 

■ naiH] 

BflaHa 

taiHaa 

■ auafli^ 1 


«d9Mf 

i44i4Q^D-l 

■i^aD^v 


+im9D 

■ KBIH] 

■(iaHa< 

‘ +ai»iMa' 

B4D<iat 1 

BaaiH 

pfltaga 


_ tpflgpg - i-mPM lABana zOUIIA LUUO lUHA. 

■DQ4«i> iiMiiaH! ^tEuta ihiaDD 




Tjnm? 


JiMfUi 






MlyHUIR 


iiimpil 


ppciLiB [ t^in flffpt# Phj n« ■ ^1, i-gji.^ E-r f r>H , p njlfr,, h Hwii-i w fh ; 

IHMTJFIC4TI4N 


NmIRf^ R aNlp Y iTfl^tri^^ltaJ tiAtIC Qwtfw 


^^»UHaB_. 


.HAQ49P 


ifliMaM 

I4MDEMK 


q4«BDW 


BlWHiaDQ 

jAAAAj3lQA_. 






■ gcBMiaB 
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«OOO0CSQO 


«OOOOGOO 


SEGMENT symmetry €H€cK 


«OQOOOQO tOOQOOOO 


<9 9000000 


90800000 


« 0000800 


oOOOOOOO 
9 008 00 bo 


90080088 


TTWOOOSW” 


*9loOOOoO‘'^01"' 




TOooOooo 


90000000 


fOOOOOOO 


rmoooootoi 


9 ««<}7656'>Os 


^008800 


eOQOOQOO 


SlMooF 


«iooo&os^oi 


«ioooooo^oi 


«t080000'^0l 


«IOOOOOQtQ| 




« 0000080 


90000008 


90000000 


VIOOOOOO'^OT ' rioooooo^oi ' »I000008'^0I 9lQ0003tT0i viOOOOOO'^Or 90000000 90000000 

IToooooiPol MooWoo^ol TIoooIiT^ol rToooooa^oi 9ioooooo^oT •Toooooo^oi iTl8l53P*o5 ?ooooooo 


sl00OOOO4>or' 


«rooods&<90i 


ftl0Qfi0Q0<90l 


ft 0000000 

PftMeWT IOaO "RTf RrCEs " 



<3e&0*5««0t 

»noonnon 




«»720&tl0»0^ 

flQQOQQQQ 




«8S900^7<»01 

Mfioaonao 


“ 


•27iSt00«03 




RZEROU$ » ^•OOOOOO^'Qa 

RZ£HQC^) » 29000000*02 















* 

/ TME»e AWe t SgcW£W Ts ANO 0 KiNcMATfC LfNKE •|TmIN THIS R£6|0n 





SE6HENT“MVH»E:R I " sewent - - 

- OTi^tf : TJTrr ^TEF "iDEETA ^ ;: 

• ITBOOCrO^Ol' riO0a0OU^Cl3 •Em0OU«*O2 tOOOOOOo 

EEOHETRT input varIabees 

. • POMOQP ^ *00 00000 •APOQOPO : 

tsOT aeuh sing thic ~ Noth Hardening Ear • eine ^wuNbEr TANJtE CpEgHNs • 2 


HA TER I AE PROPERTY TABEE 


•00000 

•00000 

• 00000 

•00000 

•00000 

•OOOOO . 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO 

•lo&oo^oa 

*00000 

•00000 

•00000 

•00000 

•00000 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO 

J^OOO^OQ-- . -- 

*00000 

_ * 00000 

•QOQOO 

ftOOOOO ~ 

. tOOOOO - — 

•OOOOO 

tiuioao. 

. *00000 

*00000-. 

•00000 

•00000 

•00000 

•00000 

•00000 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO 

*00000 

•AAAGp^pA 

• OOOOO- 

•.OOOOO— - . 

. - •OOimO_ •nnoiin . . 

^ooooo 

■OOQqO 

...^QOOOO 

-^•.00000 . . 

.. ft fiOOOO 

•aaaio^oi 

•00000 

•00000 

•00000 

*00000 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO 

AiiOOOOtflR 

-•lOOOOO 

JIOOOOO 

-- •00000 — — 

— ■*.00000^, 

-fOOOOa 

• OOOOO ^ ^ 

ft-OOOOO- 

^00000- 

•OOOOO— 


' ta©ce“oroek rni 

. -•.iiL7§.ooii?O0-..- : 

• i2Gaooo>oo •I2a3ooo'^oo 


PROBEEH i TABEE OROER JR„S :¥$, 0^ Cf.jHETA, f PH^f f 

LOaT iDENUFiCATiON CEUES^^ 

“‘“I i^^ocmb^o 4 ~'"'t&oddooo^oH ‘ ", “ “ " '/ ' 


HatR U X and Y <tNRNSpO$EP) HAGIC OUTfUT 

/doooooo •boodooo •uooQoob «booooQO «93Hbo2^bi *0000000 *0000000 tOoooooo 

^^jQOOOOO^ ^ •3.3-301ilJt0ft „^.„ji0000jtt0fi ^..^tiiOOOoaO-. , ^toaooooo ft3^lt2LlAOjtOl tOOOOOOO .-«.OflQDftflO. . 

•0000000 *0000000 *0000000 *0000000 *0000000 *0000000 *tooooootoi «ooooooo 

*0 0 0 000 0 . ft.0000000 - . ^000000 .. tooooobo- ^JJOOOOOO *:l2fl404ftt00 *22i2i40i0i._ 

• U3BS21<*01 *0000000 *0000000 *0000000 *|6292n2*0A *0000000 *0000000 *0000000 

^<200000- ~ _ *AA«X9AZ1lOO.. aOOOOOOO 1 tOoOOOOO- : *0000000. - „.tOOOOOOO.- «oaooooo 

•0000000 *0000000 •i0A6A&6400 *3307772*01 «00000o0 *0000000 *«lpD0AN7*0A *• IS7 aB92«0B , 

_.__*OOOOOdO . *0000000. «0000000 ..ftOOOOOOO — SAa9l472*>^00 »*.GNAAfi705^4.C-^ 

• OOOOOOQ *0000000 •**92AB33<I«01 «»*A9l 9H2Pt00 *OOOQOOO *0000000 *7a223fN^OA « 1 7a<t09A*0<l 



TWCi 




roHC t 


•itTirrr 


•00090011- 


•OOOOOOO ' 


~jxmtr 


MOMCNT 


"•OOOOOOO 

fOOODOOO 


~ sTtrmcss cocff icicNfs 

oftT« ri —“~“OCtT« ifr TWnri octTO t« 

~“%400^13>OV 

• 3t 1 1 4 1 s^07 •OOOOOOd 0400^ r ^0* •0000000 

— iOOOOOOO i 4041171-4407 9 0000000 ' iOOOOOOO 

'«rit07^0ffP400 4 0 000000 — VlilOorOO'mS tooooooo 

— • - 9 000 0000 40 0 000 0 0 •73^40I400 

p 9 34rtor4V0f — ^ 90 00 0 0 00 iTtO0O^1ini|4 9 ooooooo 

~»o0oooQo *rg205ot7mr7 — iinmoooor — •ooooooo 


•93H7il4340« 


9 ooooooo 


OClTO 70 

" VOOOOOOO 

••3fll4li^0T 

-90000000 

-• r370f n^04 “ 

•oooooor 

*3nm4^or 
•ooooooo 

"~V3r'»7O4S240» 


* 0CLT4 02 

“^•ooooooo 

•ooooooo 

•9*0^744S407 

•ooooooo 

-- VOOOQOQO 

““VOOOCOOO 
- 9t20Wl400 
“■90000000 " 


S£«N£NT SYNHCTRY £H7CK 


•44403§7«07 


•OOOOOOO 


•ooooooo 


•OOOOOOO 


Tnnnnnnr^'^ 


"•looooootor 


♦4i4 47 iy ^ or 


"•ooooooo 


•OOOOOOO 

"iOOOOOOo 

tooooooor 


•ooooooo 
•OoOOOOO" 
• OOOOOOO" 


•ooooooo 
— iOOOOOOO" 
•ooooooo ' 


• I 000000 YOi ^~TTI40003*^01“ — •1000000*01 


•14ygi 20*04 



• • 1 00 0 000 ^1 nooootT^oi #t 000000*0 r ~ 


nrtimao20*of 


•0000009 
VY32744l40B - 
— •-lOOOOOO^OT" 

— rtoooooom 
9igaooob*ot 


90000000 ~~ 

- “iUOOOOOO” 

--T3trW4*o^”~ 


•ooooooo 
vooooooo - 

•OOOOOOO' 


— rroooooo’*oT“" — moo2i r*o« 

9 1 000017*0 r *~T»t 000000*01 


— ^ S E 4HgH y p y^O hxtRTCEs 

90000000 

• 3H0472401 

tflOllO QQ O — : 

••4040342*OS 

- fOOQftflOO! 

•7773»28*0I 

9.0000000 ^ 

•2aot4at^oo 

iiacOoiir « I^WOooiOa o^R<riyT'.« uaTSoOoioI 


to£Y 4 a 
90000000 — 
• •3^7«2U*04 
90000000- 
•Y444tO3*0^ 

• OOOOOOO" 

93^742f2TOfc 

•OOOOOOO 

9«l4^$7i070S 

•ooooboiT 
• 0000000 """“ 
•ooooooo 

• O000000~™ 

9000D000 

•ooooooo 
•ooooooo 

94004703400 



RE«iQN NUNfER 3 

THEKE ARE t SecHEMTS *H P 0 KiNfMATtC EjNKS *|ThIN THI$ ReSION 



-- SWHtKT-miMWII— I ^ - 

BT*o Birr *TEP ; uetT* 



.au300D»BB •lOOOOOOoOO 

,8226BP0-BI •OaODOtn 





OEOHCTRY input variables 





•0000000 

*OOO0OOO *0000000 



: 







»TW6 “ThTC“ 

I4wf H H*RM*^N4 Ww law w A|wt nuno&n wr i >ui,unNa « « 

HAT&RIAU property TAflUe USEi) 

«O0OOO 

•OOOOO «00000 «00000 

•OOOOO *00000 

•ooooo 

*00000 «00000 

•ooooo 

UQSOO'tOO 

»a3obdtoo 

• OOOOO « OOOOO' « OOOOO 

•OOOOO oOOOOO *00000 

•OOOOO *00000 

*00000 *00000 

•ooooo 
_ *00000. _ 

•OOOOO *00000 

•ooooo *00000 

•OOOOO 

•OOOOO 

•OOOOO 

•OOOOO «00000 *00000 

•OOOOO *00000 *00000 

•ooooo *00000 

•ooooo *00000 

•ooooo 

•ooooo 

•ooooo *00000 

•ooooo *00000 

•ooooo 

•OOOOO 

»4ftQ0ailt04 

•OOOOO •OOOOO •OOOOO 

•aooao •oOooo •ooooo 

•ooooo *00000 

•ooooo *00000 

•ooooo 

fOOOOO 

•ooooo «ooooo 

• OQOOO.. „ •OflOOGL 

•ooooo 

. 


«l*76000t00 

TAlT-rTJRoEK PH4 OR !> V>* cR»»SfcCTlON PROPERTIES 

•AIOOOOO^OO 



»t2R3000^00 

•1283000*00 






matrix X and T (transposed) HA6IC 

Output 


' •OOOOOOO ' 

_ •opooooo 

•BOQOOOO 

•0000000 

•0000000 *0000000 

•32&3i0^t07 *0000000 

•OOOOOOO *32B333B^0I 

•OOOOOOO *0000000 

•OOOOOOO *0000000 

• *129<I3B0^0I *0000000 

•0000000 

•OOOOOOQ 

•0000000 «bbooooo 

•ooOoooo tOOOOOOO 

•OOOOOOO *0000000 

•*^4i62SBBtO^ *0000000 

•OOOOOOO 
•OO00ft00„ 

• I000000*»0l 
. •1i®|i30t00 

•OOOOOOO 

«I2-?1»6040» 

• 9i|^9?0R»Oi 

•OOOOOOO *0000000 

•A94A307«00 *0000000 

•OOOOOOO •8m823a2*08 

•OOOOOOO *0000000. .. 

•OOOOOOO *0000000 

• 1 »26Bo0^O6 _ *0000000 

•OOOOOOO 

. •oo.oofiga ... 

•OQOOOOp 
•OOOOOOQ • 

•OOOOOOO *3073770^00 

•OOOOOOO *0000000 

•1788290^00 *0000000 

♦*»*8307t00 *0000000 

•OOOOOOO 
^DOAQiLOS. 

»*32^3^03«0& 
. **304|B92*0«t 

■ • 20 e<t 828 * 0 ‘t 
•••)38ll3e^;e3 

•0000000 

•OOOOOOO *0000000 

•OOOOOOO *0000000 

•OOOOOOO 

•OOOOOOO 

•obboooo 





“ DEXT *^ lTCtTA““n OEXtft rmTh~l MXTS T 2 ■ TCCtjFi :* 6 gtT * W 2 “ THElT ^ 


^ “ -'^TORCE ■ 

—97mr&^is^urr~ 

*00011000" 

*01000000 


-♦218071 WOT — 

^irowowj 



TqrCE 

~~«ooooopo~ 

-•8627iaS^« 

“^OOoOOo 


, *60000017 

« • 8ir2TT0W05“ 

•OOOOOOO 

«*j7i8JXl0f0r“ 

rORCE“ 

“♦0000000 

♦0000000 ^ 

l|y|tiU 

•0O0O0O0 

•0000000 

•0000000 


T000000O — 

moment 

• OOOOOO'O 


•OOOOOOo 

5T50W&Wo^ 

•OOOOOOO 

•TIf0iTP0I 

•OOOOOOO 

•lioliot^of 

“'“'roHCX"";' 

■««2TSll5A‘«'i7 ■ 

•oooOdoo 

■ • 6000000 

« 6000000 


•oooOooo 

“^O0OO0O~ 

•0000000 

‘ji“ata^ n M T‘*“' 

fITrce 

• OOCTOOOO 

“w f bxstt o^IOB'” 

« 6000000 


•OOOOOOO 

•882710^^08 

• OOOOOOO. 


foscr"^ 

^'OOOO'OO 

•0000000 


•OOOOOOO 

•OOOOOOO 

•OOOOOOO 

TtWf^T^W 

•0000000 

moment 

9 0000000 


•OOOOOOo 

•25007^X«65 

•OOOOOOO 

•TW2W¥*08 

♦0000900 






iWMBiiil 

chc^e 






•7029S3S*07 

•OOOOOOO 

•OOOOOOO 


•OOOOOOO 

•OOOOOOO 

♦OOOOOOO 

♦OOOOOOO 


•loooooo^or 

iniirios^oiT 

•OOOOOOO 

•OOOOOOO 

•OOOOOOO 

•OOOOOOO 

.Oodoood 

•0600066 


“•IO0000CP01 

“•“lioWao^or 

• O2T722O807 

•6666000 

•OooOooo 

«6o666~00 

♦0O00O0O 

•0000006 


-^TTOffBoornn 

— rrsim^rrmr^ 

•T000O0CP01 

.4tt»»*V*+dS 

i066666o 

«(idoocrnr 

•0000660 

•OOOOOOO 


•1000268861 

• lottOAoona " 

•XOOOOoo^ol 

•lOQOOOO^Ol 

*0M2STf80X 

•OooOooo 

•0090000 

•OOOOOOO 


• J 660600801 

“ vtoisfSfiros^di 

• 100006680 1 

5Too06¥T8oi 

•loooooi^^^ 

•B427I6¥8pr^ 

•OOOOOOO 

•OOOOOOO 


♦ 4660060iH 

riTolOOOOVol 

•■' 7Toi®02X?oT 

• iooo66o«oi 

TToOOOOoioI 

*roo6666tol 

♦ IM8l6/8t0r“ 

,doiioSS 


•ftOOOOOO^OX 

iloo©60786l 

*roooood^t 

XTooOTXWoi 

•Toooooo^oX 

•TooOOiT^ol 

rtooooSS^di 

• 8088oO8j8ol 

' 


SE6HENT LOA& MATRICES 
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_ T phi theta 

oheoa theta Omega phI 


R ZERO base thickness 

6a¥Ha 'PHI THEU K phi 
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APPEKDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (Sl) was adopted by the Eleventh 
General Conference on Weights and Measures in i960. Conversion factors for 
the units used in this report are given in the following table: 


Physical Quantity 

U.S. Customary 
Unit 

Conversion 
factor (*) 

SI Unit 
(**) 

Length 

in. 

. 0.025h 

meters (m) 

Stress modulus 

ksi 

6.895 X 10^ 

newtons /meter^ 
(N/m^) 

Stress resultant 

Ibf/in, 

175.1 

newtons /meter 
(N/m) 

Temperature change 

°F 

5/9 

Kelvin (k) 


* Multiply value given in U.S. Customary Unit by conversion factor to 
obtain equivalent value in SI Units, 


** Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

giga (G) 

10^ 

mega (M) 

10^ 

kilo (k) 

10^ 

deci (d) 

10"" 

centi (c) 

10"^ 

mini (m) 

10”^ 
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